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ABSTRACT 


This  paper  is  concerned  with  the  Cauchy  problem  for 
linear  hyperbolic  equations  or  systems  with  variable 
coefficients.   The  main  problems  considered  are  the 
propagation  of  singularities,  focussing  of  singularities, 
and  the  construction  of  fundamental  solutions.   The 
methods  which  are  employed  correspond  closely  to  the 
methods  of  geometrical  optics  and  related  fields* 
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EXACT  AND  ASYMPTOTIC  SOLUTIONS  OF  THE  CAUCHY  PROBLEM 

lo   Introduction 

A,   Statement  of  the  problem. 

This  paper  is  concerned  with  the  Cauchy  problem  for 
linear  hyperbolic  systems  of   k  partial  differential 
equations,  each  of  order  m,   in  k  functions  (which  we 
write  as  the  vector   u(x))«   All  the  coefficients  of  the 
equations  are  generally  assumed  to  be  smooth   (C   ) 
functions  of  n+1   independent  variables  (which  we  denote 
by  the  vector   x)<,   The  chaj^acteristic  directions  of  the 
system  need  not  be  distinct,  but  their  multiplicity  must 
remain  constant  within  the  region  under  consideration* 
Our  main  objective  is  the  study  of  the  properties  of  the 
solution  in  cases  in  which  the  Cauchy  data  is  singular  in 
the  sense  of  distributions. 

We  shall  give  detailed  consideration  to  two  problems 
of  this  type  which  are  of  particular  interest  for  mathema- 
tical and  physical  applications s   that  in  which  the  data 
is  singular  along  an   (n-1 ) -dimensional  manifold  in  the 
initial  surface,  and  that  of  the  fundamental  solution  (in 
which  case  the  data  is  singular  at  a  single  point  on  the 
initial  siorfacelo   Although  the  existence  and  uniqueness 
of  solutions  with  distributions  as  Cauchy  data  has  already 


been  proved  (see  J.  Leray  [8]  and  R<.  Courant  and  P.  D.  Lax 
[ij.J),  we  shall  incidentally  obtain  a  new  proof  in  these 
two  cases  in  the  course  of  our  construction  of  the  singular 
part  of  the  solution. 

This  work  was  inspired  by  the  paper  by  R.  Courant  and 
P.  D.  Lax  [I4.J  which  deals  with  solutions  with  jump  disconti- 
nuities in  their  Cauchy  data,  and  the  fundamental  solution,— 
This  paper  extends  the  methods  and  the  results  of  Courant 
smd  Lax  to  cover  a  much  larger  class  of  problems.   Their 
methods  are  clarified  and  simplified  by  a  systematic  use 
of  the  concepts  and  results  of  the  theory  of  distributionso 

In  the  preparation  of  this  paper,  the  author  has 
benefitted  from  the  instruction,  assistance  and  advice  of 
Richard  Courant,  Joseph  B.  Keller  and  Peter  D.  Lax. 
Peter  Ungar*  made  many  suggestions  to  improve  the  manuscript. 

Bo   Method  and  results. 

The  method  of  this  paper  may  be  described  as  an 
application  of  the  methods  of  geometrical  optics  to  hyper- 
bolic initial  value  problems.   The  transition  from  wave 
optics  (represented  by  the  partial  differential  equation) 
to  geometrical  optics  (which  involves  only  ordinary 

—^Methods  of  Mathematical  Physics,  II,  by  Ro  Courant  and 
D,  Hilbert,  to  appear  shortly,  will  contain  a  compre- 
hensive treatment  of  the  problems  and  methods  of  [ij.], 
and  a  detailed  treatment  of  certain  topics  discussed  in 
this  paper.   The  propagation  of  jump  discontinuities  is 
also  treated  in  R.  Lewis  [9 Jo 
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differential  equations)  Is  effected  by  means  of  an  asymptotic 
expansion^   For  example^  we  may  look  for  formal  solutions  of 
the  wave  equation  of  the  form 

/-,  -,  \  /.   ^   _   lcJjZ^(t,x)  v^  a-^(t,x) 

(1«1)  v(t^x)   =  e   ^   ■'   2 ~ —  '  ^ 

j=0   ( ico)  J 

Here  we  neither  require  nor  expect  that  the  series  shall 
converge,  but  only  that   v   shall  be  an  approximate  solution 
for  sufficiently  large  values  of    Wo   In  other  words, 
the  partial  suin 

(1,2)  vJ{t,x)   =  e^-^^^^'^^  ^^i^   , 

j=0   (loJ)«^ 

shall  differ  from  some  solution  of  the  wave  equation  by  a 
term  of  order  less  than   CO    o 

These  conditions,  for   J  =  0,1,=,.,   immediately  imply 
relations  which  0(t^x)   and   a-^lt^x)   (j  -•  0,1 5,0,0)   must 
satisfy,,   In  particular,   |^(t,x)   must  be  a  solution  of  the 
characteristic  equation 


do) 


The  characteristic  equation  can  be  solved  by  the  method  of 
characteristics?   certain  curves  (called  bicharacteristics 
or  rays)  are  introduced,  and  0     is  determined  by  means  of 
ordinary  differential  equations  along  the  rays.   Similarly^ 
the  coefficients   a''(t5x)  may  be  determined  by  means  of 


2 

v=l     x 
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ordinary  differential  equations  (called  transport 
equations)  along  the  rays. 

Again,  it  should  be  emphasized  that  these  asymptotic 
solutions  are  generally  not  exact  solutions  of  the  partial 
differential  equationo   However,  asymptotic  solutions  are 
usually  easier  to  compute  and  manipulate  than  the  corres- 
ponding exact  solutions o   Frequently  one  or  two  terms  of 
the  asymptotic  expansion  provide  a  useful  and  convenient 
appr oximati  on. 

The  scope  of  the  methods  of  geometrical  optics  can  be 
extended  in  two  waySo   First,  such  asymptotic  expansions 
can  be  constructed  for  any  partial  differential  equation 
or  system  with  real  characteristics,  which  satisfy  certain 
algebraic  conditionSo   For  a  wide  class  of  hyperbolic  equa- 
tions or  systems,  the  Cauchy  problem,  has  an  asjmiptotic 
solution  provided  that  the  Cauchy  data  has  an  asymptotic 
expansion. 

Secondly,  we  observe  that,  in  deriving  the  equations 
which  determine   ^(t,x)   and   a-^Ctpx),   only  a  single 
property  of  the  exponential  function  was  useds   namely 


(1  L)  ^    2^  -    ^'"^ 


It  follows  that  if  we  have  a  set  of  functions 
which  satisfy  the  relations 
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(1.5)  ^f'j^^^     =    ^j«i^-^) 


then   the   expression 

00 

(lo6)  u(t,x)      =     ;^f.(0(t,x)}    a-^Ctjx)      , 

will   also  formally   satisfy  the    differential   eqiiation. 
Now   such  a   set   of   functions   is   essentially   determined  by 


o  ■ '  ■  * 


the  fiinction  f^(0)|   we  have 
(1.7) 


Moreover,  the  function  f  {0}     may  be  chosen  arbitrarilyo 

Thus,  once  0itgx)      and   a^(t,x)   are  determined,  we  have 

a  family  of  formal  solutions  which  depends  upon  an 

arbitrary  function.   Such  a  family  will  be  called  a 

generalized  progressing  waveo-^  The  function  0(tj,x)   is 

called  the  phase  function,  the  characteristic  surfaces 

0  =   consto   are  called  phase  surfaces,  the  function  ^^^0^ 

is  called  the  wave  form,  and  the  factors   a''(t,x)   are 

called  the  distortion  or  attenuation  factors., 

^  The  concept  of  a  progressing  wave  appears  in  Co-urant-Hilbert 
[Ijo  Po  Go  Friedlander  [5J  observed  that  the  infinite  series 
(lo6)  is  always  a  formal  solutiono   Friedlander  [6 J  makes 
extensive  use  of  generalized  progressing  waves  in  treating 
the  wave  equationo 
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Since  such  progressing  wave  expansions  will  not  neces- 
sarily converge,  we  must  find  an  interpretation  of  these 
expansions  which  is  independent  of  their  convergence o   We 
recall  that  the  geometrical  optics  expansion  (lol)  can  be 
Interpreted  as  a  limit  of  oscillatory  solutions  as  the 
frequency  approaches  infinity o   The  behavior  of  a  function 
at  infinity  is  closely  connected  with  the  smoothness  of  its 
Fourier  transformo   Thus  there  is  a  link  between  higlily 
oscillatory  solutions  £ind  rapidly  varying  or  discontinuous 
solutions.   If  the  wave  form  f  (jZl)   has  a  discontinuity 
in  one  of  its  derivatives,  we  observe  from  the  defining 
relation  (lo5)  that  f A^)      becomes  arbitrarily  smooth  if 
j   is  chosen  large  enough.   Thus  if  we  set 

(1.8)  u'^(t,x)  =     m  f,(i^(t,x))  aJ(t,x)   , 

then  u  (tpx)   differs  from  an  exact  solution  by  a 
function  which  becomes  arbitrarily  smooth  (not  arbitrarily 
small)  as   J  —-*>  oo  o   This  interpretation  is  valid  if  we 
choose  an  arbitrary  distribution  for  the  wave  form   f 
Thus  the  primary  usefulness  of  progressing  wave  solutions 
is  the  following  I   by  means  of  a  finite  number  of  terms  in 
the  expansion,  we  can  transform  a  problem  with  singular 
Cauchy  data  into  a  problem  with  arbitrarily  smooth  Cauchy 
data.   It  is  noteworthy  that  such  progressing  wave 
expansions  are  determined  by  means  of  ordinary  differential 
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eqxiations  along  the  rayso 

Because  of  the  great  formal  similarity  between 
generalized  progressing  waves  and  asymptotic  expansions, 
it  will  be  convenient  and  not  misleading  to  call  general- 
ized progressing  waves  asymptotic  solutionse 

The  phenomenon  of  caustics  furnishes  an  example  of 
the  fruitfulness  of  the  analogy  between  asymptotic  expan- 
sions of  geometrical  optics  and  generalized  progressing 
wave  So   Caustics  appear  where  the  rays  of  geometrical 

optics  form  an  envelopeo   At  a  caustic  there  is  a  change 

3/ 
in  the  asymptotic  behavior  of  the  solutiono—  Thus 

instead  of  a  leading  term  of  the  form 

(lo9)  V  ^   e^^^  a(t,x)   , 

we  may  have 

(lolO)  V    ^uf   e^^^  a(t,x)   , 

where   a  is  a  positive  rational  numbero   In  physical  terms, 

the  intensity  of  the  field  becomes  infinite  at  the  caustic 

as  the  frequency  increases.   There  may  also  be  a  shift  in 

phase  after  a  ray  passes  through  the  caustic o   These 

phenomena  of  geometrical  optics  have  their  counterparts  in 

the  behavior  of  generalized  progressing  waves «   We  shall 

see  that,  at  a  caustic ,  the  wave  form  is  al teredo   Instead 

of  a  leading  term  of  the  form 
3/See  Ic  Kay  and  J,  Be  Keller  L9Jo 
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(1.11)  u  --  f^(0)  a(t,x)   , 

we  may  have 

(lol2)  u  -v  ^_J<0)    a(t,x)   , 

where   f    is  an  a-th  derivative  of  f   ,   and   a  is  a 
positive  rational  niomber.   The  shift  in  phase  which  occurs 
in  optics  corresponds  to  an  alteration  in  the  wave  form 
after  passing  through  a  caustic. 

Although  all  of  the  above  considerations  are  valid 
regardless  of  whether  or  not  the  formal  expansion  converges, 
we  shall  see  that,  provided  that  certain  analyticity 
requirements  are  satisfied,  the  formal  expansions  do 
converge.   In  particular,  we  shall  see  that  the  Riemann 
function  (a  fundamental  solution)  can  be  constructed  by 
means  of  convergent  progressing  wave  expansions,  for  linear 
hyperbolic  systems  with  analytic  coefficients.   In  a  sense, 
this  result  is  an  extension  of  the  method  of  Hadamard  to 
more  general  hyperbolic  equations.   Since  progressing  waves 
can  be  constructed  by  means  of  ordinary  differential  equa- 
tions and  algebraic  processes,  hyperbolic  equations  can 
thus  be  solved  by  means  of  a  recursive  system  of  ordinary 
differential  equations. 
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C.      Examples. 

The  hyperbolic  equation  easiest   to   deal   vrith  is   the 
wave    equation  with  two    independent   variableso      All   solutions 
of   this  equation  can  be   written  in  the   form 

(1.13)  u(x,t)      =     P(x+t)    +   G(x-t)      , 

where   F  and   G  are  arbitrary  functions  of  a  single 
variable.   Such  solutions  are  called  progressing  waves. 
If  we  prescribe  Cauchy  data 

u(x,0)   =  f(x) 

d.li].) 

u^(x,0)   =  g»(x)   , 

then  we  can  immediately  express  the  solution  as 

(1.15)    u(x,t)  =  i  [f (x+t)+g(x+t)+f(x-t)-g(x-t)j  . 

If  we    do   not  require   that      f      ajid      g     be   dif ferentiable 
in  the   usual    sense,    but    only  that   they  be    distributions, 
exactly   the    same   formulas   give    the   solution  for   singular 
initial   data.      If     v(x)      is   any     C*^    function  with  compact 
support   then  it   follows   immediately  that 
.00         ^00 


u(x,t)    (v      -v.. )    dx  dt     =     0 
'  XX      tt 


-00       'J  -00 


Thus   u(x,t)   is  a  "weak  solution"  of  the  wave  equationo 
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However,  It  Is  unnecessary  to  Justify  our  procedure  by  means 
of  such  integrals.   On  the  contrary,  if   f  and  g  can  be 
represented  as  limits  of  twice  continuously  dif ferentiable 
functions  in  some  topology  which  respects  addition  and 
multiplication  by  smooth  functions,  then  our  formula   (l.l5) 
shows  that   u   is  a  limit  of  twice  dif ferentiable  functions 
in  exactly  the  same  sense.   This  interpretation  of  weak 
solutions  is  more  direct  than  the  usual  one,  and  more 
meaningful  for  physical  applications.   An  important  result 
of  this  paper  is  that  this  direct  interpretation  can  be 
extended  to  hyperbolic  systems  of  any  order,  and  with  any 
n\jmber  of  independent  variable  So 

Another  example  of  a  progressing  wave  solution  is 
furnished  by  spherically  symmetric  solutions  of  the  wave 
equation  with  four   independent  variables.   These  may  be 
expressed  in  the  form 

(1.16)  u(r,t)      =     F(t^r)-F(t-r) 

where  again  P  is  an  arbitrary  function.  Spherically 
sjHTmietric  solutions  of  the  wave  equation  with  any  even 
number  of   independent   variables   have   the    foiW^^ 

(n-lj.)/2      .  , 

(1.17)  u(r,t)      =     7"^     hJ(r)    [(-l)JF,(t+r)-P,(t-r)J    . 

j=0  J  J 


S^ 


See   R.    Courant   and   D,    Hilbert    [1],    p.    ipL5. 
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Here  P  (s)   is  an  arbitrary  function, 

^P^.(S)     =    Pj-_1(S)   ,  (j   =   1,    coo)   , 

and  the  functions  h'^(r)   are  regular  for   r  >  0  o   These 
solutions  are  called  "progressing  waves  of  higher  degree" 
(hoherer  stufe)o   Corresponding  to  these  solutions,  weak 
solutions  can  be  defined  and  given  the  same  direct  inter- 
pretation as  in  the  case  of  two  independent  variables. 

A  noteworthy  fact  about  the  solution  (lol6)  is  that 
the  differentiability  properties  of  the  Cauchy  data  are 
preserved  in  the  solution  except  along  the  line   r  =  0  » 
This  line  is  the  locus  of  the  singularities  of  the  charac- 
teristic surfaces  r+t  =  consto   Such  loci  of  singularities 
are  called  "causticSo"  We  shall  see  that,  for  any  hyper- 
bolic equation  in  any  number  of  independent  variables, 
differentiability  properties  of  the  initial  data  are 
preserved  in  the  solution  except  at  caustics  of  certain 
characteristic  surfaces. 
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2.   Generalized  Propiresslng  Wave  Solutions  of  Systems 
of  First  Order  Sq\xations« 

In  order  to  Illustrate  the  method  to  be  followed  In 
the  more  general  case  of  systems  of  higher  order  with 
multiple  characteristics,  we  shall  first  construct  a  family 
of  formal  solutions  of  a  system  of  first  order  partial 
differential  equations  with  distinct  characteristics.   Let 
X  =  (x  ,x  ,o..,x  ),   and   let  u(x)   be  the  solution  vector 
with  k  components.   Consider  the  system  of   k  linear 
first- order  equations: 


(2.1)  <7^(u)   =  A^  ^  + 


Bu  =   0 


Here  the   k  x  k  matrices  A   and  B   are  infinitely 

differentiable,  and  the  index   v   is  to  be  summed  over 

from   0   to   n  whenever  it  is  repeated  in  a  single  term. 

Provided  that  the  planes  x  =  const.   are  nowhere  tangent 

to  a  characteristic,  we  can  bring  the  system  (2.1)  into  a 

form  where  the  matrix  A   is  the  Identity.   We  assume 

that  this  has  been  carried  out. 

We  shall  construct  a  formal  solution  of  (2.1)  which 

depends  upon  an  arbitrary  function  of  a  single  variable: 

f  (s)  .   First  we  define  the  functions  f ,  ,f -,,...  ,f  .... . 
o  1'    £i'         '    y 

in  such  a  way  as  to  satisfy  the  conditions 

^  f.(s)   =   fj.i(s)  (j  =  1,2,...). 
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Tlien  we   assume    that   the    solution  has   the    form 

00  . 

(2.2)  u(x)      =     m  ^^    (0(x))    &Hx)    . 

f   is  called  the  wave  form,   ^(x)   is  the  phase  function, 
and  a''(x)  are   the  distortion  coefficients.   Here  the 
coefficients   a*^   and  the  function  0     are  to  be  deter- 
mined.  Next  we  substitute  the  expansion  (2.2)  into  the 
differential  equation  (2.1) s 

(2.3)  ^u  =  f  ,(0)[AXa°]  +  21:  f.(^)[AV  a^*"'^+o^(a^)]  =  0, 

In  (2.3)  and  henceforth  we  employ  the  notations  0  =  S0/dx    , 

In  order  to  determine  the  coefficients   a"^(x),   we 
require  that  the  coefficient  of  f .   vanish  identically 


for  each  value  of  2*        Thus  we  have 
(2.k)  a''  0^  a°   =   0  , 

(2.5)  a''  0^   aJ""^  +  ^{&h      =0      (j  =  0,1,...) 
Hence,  from  (2.14.), 

(2.6)  det  (A^'^y)   =  Ai0^)   =  0,    ^ 

^  The  calculations  in  the  remainder  of  this  section  are 
similar  to  those  in  R.  Courant  and  P.  D.  Lax  \.\\.\ ,  and 
identical  to  those  in  P.  D.  Lax  [10 j. 

— '  A  vector  subscript  will  always  denote  the  gradient 
with  respect  to  that  vector. 
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Equation  (2«6)  Is  the  characteristic  equation  corresponding 
to  the  system  (2.1).   It  follows  that  the  phase  surfaces 
0ix)   =   const,  must  form  a  family  of  characteristic  surfaces 
of  equation  (2.1). 

A  second  consequence  of  (2.[|.)  is  that   a   must  be  a 
right  nullvector  of  the  characteristic  matrix  A  =  A  0'  « 
For  convenience,  we  assume  that  A  has  rank  k-1.   This 
restriction  will  be  removed  in  Section  !+<,   Then  we  may 
define  the  right  nullvector  R(x)   by  the  conditions 

(2.7)  AR  =  0  ,      lR|  =  1  . 

Consequently     a       must  have   the   form 


(2a8)  a°  =  cr°  R 


g 


o 
where      cr         is   a   scalar   factor   still   to  be   determined. 

Upon  setting      j   =   0     in  equation    (2.5),  we   have 

(2o9)  A   a^   +     -^(a°)   =  0    . 

We  may  consider  (2.9)  as  a  system  of  linear  equations  for 
the  vector   a   ,   where  the  coefficient  matrix  A   is 
singular.   This  system  of  equations  has  a  solution  if  and 
only  if 

(2.10)  L  c^(a°)   =   0  . 

Here   L  is  a  left  nullvector  of  A  s 

LA  =  0  o 
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Substituting   (2.8)    into    (2.10)   yields   the   following 
ordinary  differential   equation  for      o"       ; 

(2.11)  lA^  R  cr  °  +  L  c^(R)   cr°     =     0   , 


We  shall  see  presently  that  the  direction  of  differenti- 
ation in  (2.11)  is  bicharacteristic.   Provided  that  (2.11) 
is  satisfied,  we  can  solve  (2.9)  for   a  ,   which  will 
then  be  determined  modulo  R.   Therefore  we  may  write   a 
in  terms  of  a  particular  solution  h   and  a  scalar   o^ 


in  the  fointi 


a        =     h     +  <r      R   . 


Now  we  proceed  by  induction  on   j  .   Let  us  assume 
that   a°,...,a''   have  been  determined  in  such  a  way  that 
the  coefficients  (2.I4.)  and  {2,S)    of  f  _^,f ^, . . .  ,f  .  ^  are 
identically  zero,  and  that   a*'   Is  determined  modulo  R: 

a''   =  h'^  +  cr  "^  R  . 

We  set  the  coefficient  of  f .   equal  to  zero  (see  eq.  Z,S)'» 

(2.12)  A  a^"^^  +  tuh      =   0  o 

The  unknowns  in  (2.12)  are  the  scalar  c  ^      and  the  vector 
h*'    (=  a*^   mod  R).   The  necessary  and  sufficient  condition 
for  the  system  (2.12)  to  have  a  solution  is 

(2.13)  L  (^Uh      =     0  . 
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Hence, 

(2,114-)  L  A^  R   crj    +  L,y!(R)    cr  ^      =      -L    tilih    . 

This    is    an  inhomogeneous   ordinary   differential  equation 
for      cr^,      which  has   the    same   left-hand   side   as    (2oll), 
with     o"  "^      in  place   of      o-  °,      Once      o-  ^      has  been  deter- 
mined  J    we    can   solve    (2.12)   for     h*^      .      Thus    the   induction 
is   complete,    and  the    coefficients      a^      can  be   determined 
successively. 

Conversely,    suppose   that   a   f^anction     ^     and  vectors 

o      1 
a   ,8.   5  o  «  o ,      can  be   determined  which  satisfy  the   relations 

(2o6)  det    (aV^)      =     0    , 

(2.8)  a°     =      a-°  R    , 

(2ol2)  A  aJ"^-"-   +     c£{a.h      =      0   ,       (j   =  0,1,...). 

Then  our  calculations  show  that   u(x),   given  by  (2.2),  is 
a  formal  solution  of  the  system  (2.1).   u(x)   is  uniquely 
determined  if  initial  values  for  )2S(x)   and   cr  (x)   are 
assigned. 

¥e  now  prove 
Lemma  2ol;   The  direction  of  differentiation  in  equations 
(2,11)  and  (2.1L|.)  is  bicharacteristic  and  ^(x)   is  constant 
along  a  bicharacteristic  ray. 

The  bicharacteristic  direction  corresponding  to  a 
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characteristic    surface      0  =  consto      is   given  by 

^     =     oMif-  (v  =  0,...,n), 

where  c(r)  is  some  constant  of  proportionality.  We  can 
express  this  ratio  in  a  more  convenient  fonn  if  we  deter- 
mine 0       as  a  function  of  0^,  ,<,  o  ,0       from  eq.  (2e6)s 

(2.15)  0Q   +  ^(0^,'",0J      =  0  . 

Note  that  the  function  '>^i0-,^  -  <> '  i,0^)      is  homogeneous  of 
degree  1  in  its  arguments,  while   ^i0   , ' ' » >0^^      is  homo- 
geneous of  degree   k  .   Equation  (2.15)  is  more  convenient 
to  work  with  than  equation  (2.6),  since  the  solution  of 

(2.15)  is  uniquely  determined  by  its  initial  values.   A 
solution  of  (2.6)  need  not  be  uniquely  determined  by  its 
initial  values  because   £>.(0^,j2!^, . .  .,0^)   is  a  polynomial 
of  degree   k   in  J^f  .   A  solution  of  (2.15)  is  uniquely 
determined  by  its  initial  values  because  0^      is  given 
explicitly  as  a  function  of   (0,  , , , .  ,j2f  ). 

Now  the  parameter  "^   may  be  defined  by  setting 

(2.16)  i^  =  1.  1^  =  #;      (,  =  l,...,n). 
Upon  differentiating  (2.7)  with  respect  to  0  ,   we  obtain 

If  we  multiply  this  equation  by  L,   we  obtain 
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(2.17)  L  A^"  R   =  i^  L  A°  R   . 


Upon  combining  this  result  with  (2.16)  we  obtain 


,  V 
(2.18)  L  A^  R  =  L  A°  R  g-     ( v  =  0,1, .  ..  ,n) 

Equations  (2ol8)  show  that  the  direction  of  differentia- 
tion in  (2«11)  and  (2. 114.)  is  bicharacterlstiCo   Furthermore 
the  homogeneity  of  the  function  X   implies  that 

Hence  0     is  constant  along  a  bi characteristic  and  the 
lemma  is  proved. 

We  may  siommarize  the  results  of  this  section  in  the 
following 

Theorem;   u(x)  ,   given  by  (2<,2),  is  a  formal  solution  of 
(2,1)  for  arbitrary  choice  of  f   ,   provided  that 

lo  the  phase  function  0   satisfies  the  characteristic 

equation. 
2,  the  scalar  factors   er  -^   satisfy  the  ordinary 
differential  equations  (2,11)  and  (2oli|.)  along 
the  bi characteristics  which  generate  the  surfaces 
0   =  consto 
3«  the  distortion  coefficients   a*^   satisfy  the  linear 
algebraic  equations  (2«12), 
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3o   Generalized  Propjressinp;  Wave  Solutions  of  the  Cauchy 
Problem  for  Systems  of  First  Order. 

In  this  section  we  shall  use  the  formal  expansion  of 
Section  2  to  construct  the  formal  solution  of  a  Cauchy 
problem  with  initial  data  involving  an  arbitrary  functiono 
We  denote  x°  by  t   and   (x  ,...,x^)   by  y  o  We 
consider  the  totally  hyperbolic  first  order  system: 


(3.1)  o^u  E  A^^^ 


^x 


Bu  =  0  , 


where  the   k  x  k  matrices  A   and  B   are   C   ,   A  =1, 

and  the  hyperplanes   t  =  const.   are  spacelike:   i.e.   all 

real  linear  combinations  of  the  matrices  A  ,..,,A   have 

k  real  and  distinct  eigenvalues.   The  more  general  case 

of  multiple  eigenvalues  will  be  treated  in  Section  5. 

Let  us  assume  that  initial  data  for  (3.1)  is  prescribed 

in  the  form  of  a  wave : 

oo  . 

(3.2)  u(0,y)   =  JZ  ^A0(j)^    b^^y)  . 


J=0   ^ 


where 


1,  f   is  an  arbitrary  function,  and 

^  fj(s)   =  fj._i(s)  ,    (j  =  0,1,...), 

2.  ^  is  a  C°°  function  of  y  ,   and  the  gradient 
of  0     with  respect  to   y  does  not  vanish. 
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3.   The  vectors   b'^(y)   are   C°°   functions  of  y. 

7/ 
First  we  shall  construct   k  distinct  functions-'-/ 

0  (x)   which  satisfy  the  characteristic  equation 

(3.3)  det  (A^i2^J)   =  A(j^J)   =  0,   (a=l,.,.,k), 

with  the  initial  condition 

(3.1^)  i2^^(0,y)  =  i2(y)  ,       (a  =  l,,..,k). 

In  geometrical  terms,    we   pass      k     families   of   characteristic 
surfaces    through   the   phase    surfaces   of   the    initial   data. 
This    is   equivalent    to    solving   the    following   set   of   ordinary 
differential    equations   along   the  rayss 

with  the   initial   conditions 

x^^(O)      =      y^^   ,      j^^(o)    =  -i-  jZ((y)       (ix=l,...,n)j 

,  c)t  7 

(3.6)^  Ma=l,..0  5kJL 

x^(0)      =      0,     ^(0^(0),    0j{O),,. .,0^(0))    =  0  ) 

The   equation 

A(0j(O),    0j(O),     ..o,    0"(O))    =    0    , 

which  appears    in    (3o6),    is   the    condition  that      -  ^    (O)      be 
an  eigenvalue    of   the    matrix      0-.(O)   A      +    „..    +0    (0)   A    , 

—'    This   is   possible    since   the   plane      t   =  0      is    spacelike. 
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Since  we  have  assumed  that  the  plane   t  =  0   Is  space - 
like,  there  are   k  distinct  roots  of  the  eigenvalue 
equation.   Hence  there  are   k  distinct  solutions  of  the 
system  i3'S)   with  initial  condition  (3<>6),  and  k 
distinct  solutions  of  equation  (3»3)  with  the  initial 
condition  (3oi|.)o   We  assume  that  there  is  a  number 
t-,  >  0   such  that  the  solutions  of  (3<>3)  are  regular 
for   0  <  t  <  t-,  o   Such  a  number  always  exists  if  the 
Cauchy  data  (3=2)  has  compact  support » 

Next  we  consider  a  formal  expansion  u(x)   consist- 
ing of   k  expansions,  each  in  the  form  of  a  generalized 
progressing  wave.   Then  we  determine  the  initial  values 
of  the  coefficients  crHx)      in  the  expansion  from  the 
initial  data  (3,2),,   We  write 

00   k      „  „. 

(307)  u(x)   =  y^  y~  fA0°-{x))   a^J(x)  , 

j=0  a=l   J 

To  each  of  the  functions  0       are  associated  unit  vectors 
R   and  L   such  that 

A^  i2f^  R^  =   0  , 

(308)  V,        a  =  1,  o  e  <,,k, 
L°^  A^  jZ^;  =   0  o 

We  define  the  scalars   o   °(x)   by  requiring  that 

(3olO)         a^'^lx)   =  er^°(x)  R^'lx)  ,     (a  =  l,o..,k) 
Thus,  from  (3.7), 


2k   - 


(3.11)  u(0,y)  =  fA0ij))['ZZ   o  ''°(0,7)R''(0,y)]  + 

°       a=l 

j=l   -^ 

Upon  identifying  the  coefficients  of  f^  in  (3.11)  and 
(3.2)  we  obtain 

(3.12)  ^   0-''°  (0,y)  R''(0,y)   =  b°(y)  . 
a=l 

Since     ji   ,^  ,o.,,j2^        are   distinct,    the    vectors 

1  k 

R   ,,.o,R        form  a   complete    set,    and   the    factors 

cT      jcr       jcoojC  can  be  uniquely  determined.      Once   the 

initial   values      cr  "^"^(Ojy)      have  been  determined,    we    can 

define    the    functions      C    ^      along  the   rays  which  start   at 

(0,y)      by   integrating  the   ordinary  differential   eqviations: 

(2.11)  lVr^(^*^°  +  L^ot(R^)    cr  '^°      =      0        (a  =  l,...,k). 

Then  (2.10)  yields   a°°(x)„   This  choice  of   a'^°(x)   insures 
the  compatibility  of  the  systems 

(2,9)  a''  j2fj  a''^  +  oZ^(a''°)   =   0,    (a  =  l,..,,k). 

Therefore  we  can  determine  the  vectors   a    within  a 
multiple  of  R^s 

a^^   =  h^l  +  cr  "^l  R^  . 


Retiirning   to   the   initial    conditions,    we    identify  the 
coefficients   of     f,      in    (3.2)    and    (3.11)    and   obtain 
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a=l 
or 

(3.13)  YZ  ^""^  R^     =     b^  -  ^h^l   . 

a=l  a=l 

In  general,    we   determine      cr'^Ho^j)      by   solving   the 
system  of   equations 

(3.11^-)  zz  ^""^  r""  =  bJ  -  5i:  h^^  , 

a=l  a=l 

Then     o^  *^*^      are   defined   along   the    rays   emanating   from 
(0,y)      by   integrating  the   equations 

(3el5)    lVrV/J   +  L^^(R^)    cr^J   =   -L^.^Ch'^J),       (a=l,.,.,k) 

h  "^     can  be  determined  by  solving  the  algebraic  system 

(3«16)         A^  0^   a^  j"^^  +  c^ia'^h   =  0   ,      (a  =  l,.c,,k). 


Thus  we  can  obtain  all  of  the  terms  of  the  series  (3.7). 
In  summary,  we  recall  that  the  method  of  Section  2 
enables  us  to  construct  generalized  progressing  wave  solu- 
tions of  any  linear  first  order  system  of  partial  differen- 
tial equations  with  one  or  more  real  distinct  characteristics, 
In  this  section  we  have  considered  a  totally  hyperbolic 
system  and  have  assigned  Cauchy  data  in  the  form  of  a 
generalized  progressing  wave.   We  have  shown  that  if  the 
solution  is  assumed  to  consist  of  a  sum  of  generalized 
progressing  waves,  then  it  has  a  imique  decomposition  into 
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a  sum  of  such  waves.  Each  of  these  waves  has  a  family  of 
characteristic  surfaces  of  the  system  for  phase  surfaces. 
These  formal  solutions  will  be  discussed  in  Section  6. 
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[|.,   Generalized  Progressing:  Wave  Solutions  of  Systems 
of  Hig:her  Order o 

In  this  section  we  shall  extend  the  method  of 
Section  2  to  cover  the  case  of  a  system  of  equations  (or  a 
single  equation)  of  higher  order,  where  the  characteristics 
need  not  be  distinct. 

In  order  to  make  ideas  clearer  and  to  save  writing, 
it  will  be  convenient  to  use  the  following  notation,  which 
is  due  to  Lo  Schwartz.   Let 


^    ^0*°°°*  V  n-^  * 
OX.  dx 

dx        c)x 

a  =    [a   , . oo .a    )    s      (a.   are  integers  >  0)  , 

a   a,       a 
a  _   o   1       n 

P   -  Pq  P]^   » » o  P^   • 

We  shall  have  to  make  use  of  the  rule  for  differentiation 
of  a  product  of  two  functions s   u  and   u  ,   This  rule  can 
be  expressed  in  a  simple  form  if  we  represent  D  as  the 


^>wl 


/s 


sum  of  the  operators  D   and  D,   where  D  acts  only  on 


■^  "^  /\ 


u,   and  D  acts  only  on  u.   If  Q(p)   is  a  polynomial  in 
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Pj,      then  algebraic  manipulation  shows   that 

Q(p+p)      =     Q(p)    +  G   (p)  •  p  +   .0.    ; 
consequently 

Q(D)(u:  u)      =      Q(D+D}    (u   u) 

=      Q(  d')   u   u   +  Q    (d')  0  D  u  u  +    o  o  o 

P 

=    (Q(D)il)u   +    {O(D)n)  o  D  u  +    ,.o    . 

P 

Nov/  let  us  consider  a  k  x  k  matrix  Q,  whose  elements 
are  polynomials  in  p  with  coefficients  which  are   C°° 
functions  of  Xo   We  write 

Q(p)   =  A(p)  +  B(p)  , 

where     A(p)      is   homogeneous   in     p      of  degree  m,      and     B(p) 
contains   only  polynomials    of   degree   less   than     m.      Let     u 

be   a   vector  function  of     x     with     k      component  a »   We    shall 

th 
be    concerned  with  the   linear     m  oz'der   system  of  partial 

differential    equations? 

(liol)  o^(u)   3  Q(D}   u  5=  A{D)u  +  B(D)u  =   0    . 

We  make  the  following  ass-umptlons  about  the  matrix  A: 
lo  If  the  characteristic  polynomial  has  a  root  of 
multiplicity  r,   then  the  characteristic  matrix 
has   r   linearly  independent  right  and  left 
nullvectorso 
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2o  A    is  proportional  to  the  identity. 

3e  The  multiplicity  of  each  root  of  the  characteristic 
equation  is  constant o   This  condition  implies  that 
p    can  be  determined  as  a  regular  function  of 
(,Xj,p-|  J,  o .  o  jP  )   from  the  characteristic  equation. 

As  before,  let   f  (s)   be  an  arbitrary  scalar  function 

and  let   f  .(s)   satisfy  the  relations 
0 


^  fj(s)   =  ^j-i^s)  (j  =  -^'•l^ 


o  o  »  /  o 


If  we  first  set 


u(x)   =  f^(^(x))  a(x)  , 

then 

\ 

c€u  =  f  „(j2f)[A(D0)J  +  f   .t(0)[A^(D0)o  Da  +  CaJ  +  ..o, 
—m  ='m~x     p 

where   C   is  some  matrix  factor.   In  order  to  make  our  later 

formulas  more  explicit  we  shall  calculate   C.   Since   C 

and  all  other  coefficients  In  the  expansion  are  independent 

of  f  5   they  can  be  determined  easily  by  a  particular 

choice  of  f  ,   If  we  set 
o 


(i;.2)  g„(s,x^)   =  i  (3=i^(x„)r 


then 


^*^'^o'  -  51 


^  Sm"^<"''^o' 


» 


3 

=   5 
m 


^=^o 


30 


Thus   C   can  be  determined  from  the  condition  that 


(l4-0) 


Ca 


x=x 


x=x 


,{g^_-L(i2^(x),x^)a(x^)) 
'o  "'  "o 

Now    if  we    assume    that   the    solution   of    (14..  1)   has 

the   form 

00 


u{x)      =     IZ  ^.    (0(x))    a^(x)    , 


J=0 


J 


then 

ik^S)       c^u     =     f_^(0)[A(D0)a°]   + 

+  f_^^^(iZJ)[A(D0)a-'-+A    (D{2r)  •  Da°   +  Ca°]   + 


00 


J+1. 


1  J+Po  J+K"J 


+  51  f^i+l+j^^^'^^^^^^^     -'+A    (D^}»DaJ+Ca-J+KJ], 
J      XI  ^  T  1+1-m        j+2-m  1-1  , 


where      K^      depends    only   on     a 
is    in  fact   given  by   the   formula 


, . . » ,  a'   , 


(1^,6)  KJ(a°,...,a^*"^) 


i-1 


4,  iz  ^,^..^/0^^u^)^h^)]\ 


'X=X. 


As  before,  we  determine  the  function  0     and  the 
vectors   a  ,a  ,.,,   by  setting  the  coefficient  of   f  .(^) 
equal  to  zero  for   j  =  -m, . • »  .   Thus  we  have 


(l|.c7) 


k{D0)    a°  =  0  , 


(1|.,8)    A(Di2f)a^'^-^  +  A  (Dj?^)  .  Da^+Ca^+K^  =  0   (j  =  0,1,...), 
Hence,  from  (I|^7) 
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(i|o9)  det    (A(D0]j      E       a(D^)      =     0    . 

• 

This  is  the  characteristic  equation  corresponding  to  the 
system  (Uol)<>  It  is  a  nonlinear  first-order  partial 
differential  equation  for  the  function  ^(x).  Note  that 
A(D0)  does  not  explicitly  involve  the  function  0,  and 
A(D0)  is  homogeneous  of  degree  mk  in  D0  o  IT  0 
satisfies  the  characteristic  equation,  then  the  si^rfaces 
0{x)    =   const o      are    a  family   of   characteristic   surfaceso 

If   the   rank  of     k{T)0)      is      k-r,      let      R^Cx)      and 
L'^lx)      (PsY  ~  lj,o.oj,r)      be   linearly  independent    unit 
right   and  left   nullvectors   of      A(Dj3)o 

Prom  equation    (L(^7)j.    we   see   that      a       must  be   a 

1  r 

linear   combination  of     R   ^„oo^R'      % 

(iiclO)  a°     =     orV  +   0,0    +   a-°R^     ^      o-°rP. 

1  r  p 

Setting   j  =  0   in  equation  (l].o8)j  we  have 
(ij.,11)         Aa""-  +  A  «  Da°  +  Ca°  =  0  , 

or 

(1^.12)         [Aa^  +  o^(g^  .i0,7i    )    a°]|     =  0. 

'x=x 
o 

This  equation  has  a  solution  if  and  only  if 

L'^^(gj^_l(i2^,x^)a°)|     =0,       (Y  =  l,.o.,r) 


'x=x 

o 


Hence 
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Now  we    show   that   otop   assi:mptlons   about    the    system 
(L|.ol,)    imply   that    (i^ol3)    is    a   system  of   ordinary  differ- 
ential  equations.^    and   the   direction  of   differentiation 

is  bicharacteristico   That  is^  if  we  set 

Po 

then   L'  A        R^      can  be  written  in  the  form 
Pv 

(i^oii^)  lY  a   rP  •=  A^^^  p'yP  , 

Py 
for  all  Y   and   p  ,   and  for   v  >  Oo   It  follows  from 
our  assumptions  1  and  2  about  system  (.I4.0I)   that   p'^^   is 
a  nonsingular  matrix » 

In  order  to  obtain  ([[.ol^Jp  we  write 

in  a  certain  neighborhood  of  a  particular  root  of  (l4.o9)» 
The  bicha.racteristic  direction  and  the  parameter   'gr  can 
be  defined  by 

li^olbj       ^—  -  1  ,   -^^  -  j^       {^i  -  l,...,n). 

X   is  homogeneous  of  degree  1  in  0,j,<.».j,0   ^   consequently 


33  - 


Hence 

and  0     l3  constant  along  bicharacteristics. 

As  in  the  proof  of  lemma  2ol,  we  differentiate  the 
equation 

A(DjZl)  rP   =  0  , 

with  respect  to  0       and  obtain 

^r.     jt^"^      =     K     R'^  +  ArP   ,      (t.  =  l,.,.,n). 
Hence 

L^R^l^  =  ^X  R^  .  (ti.  =  l,.o.,n) 

and  from  (ij..  16), 

(i4.ol7)     L^A  rP   =   (L^A   RP)  1^  (v  =  0,1, ...,n), 

Py  ^O 

for  all   p   and  y»   The  relationship  expressed  by 

equation  (ij..  1?),  which  greatly  simplifies  the  treatment 

of  multiple  characteristics^  was  first  noticed  by 

Ro  Couranto— '^ 

It   follows   from    (i4.ol3)    and    (lj.oll|)    that 

(1|.18)  pYP  ^  +  pYP   crp°     =     0  (y  =  l,...,r), 

where 


■2/   There    is    a    similar  remark  by   P.   Ungar,    quoted   in 
R,    Lewis    [9],    po    577» 
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,yP 


P^P(x^)   =  L^-^-(g^.i(0,^o^  «'^ 


x=x 
o 


Equations  ([4..  18)  form  a  system  of  ordinary  differential 
equations  for   cj:;^,  » « ,  ^cr  °^   since  each  term  is  differ- 
entiated in  the  same  direction.   If  these  equations  have 
been  satisfied,  we  can  determine   a   modulo  a  right 
nullvector  of  A  from  (.i4.oll)s 

a   =  h  +  cri  R^  , 

Now  we  proceed  by  means   of   an  induction  on      jo      Assume 

that      a   jooo<,a^      have   been   determined  in   such  a  way  that 

the    coefficients   of      f      .f      ,,„..., f,  in   (ko5)    are 

-m'  -m+i    *  J~i'i 

identically  zero,  and   a-^   is  determined  modulo  a  right 
nullvector  of  As 

(i|.  19)  a^  =  h^  +  o'^  RP  , 

We  set  the  coefficient  of  f^+i_n,   In.  il^oS)    equal  to  zero: 

(I|.o20)         Aa^"^-"-  +  A  o  Da^  +  Ca^  +  K^   =  0  . 

P 

¥e  determine   a"  "^   by  means   of  the  compatibility  equations 

for  ([|.o20)s 

dcr  J 
(i4..21)    p^P  -g|r-  +  P^^  cr^      =     -L^lA  •  Dh^°  +  Ch^  +  K^ )  , 

Once  ifr,        have  been  determined,  we  can  solve  ([|.«20)  for 
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h^   o   Thus  the  induction  is  complete,  and  the  coefficients 
a^   can  be  determined  in  succession. 

We  conclude  that   u(x)j   given  by  {l\.»l\.) ,    is  a  formal 
solution  of  (i|.ol)  for  arbitrary  choice  of  f  ,   provided 
that  the  following  conditions  are  satisfied. 

lo   The  function  0     satisfies  the  characteristic 

equation. 
2o   The  scalar  factors   cr„*'   satisfy  the  ordinary 
differential  equations  (i].ol8)  and  il\.<, 21)    along 
the  rays  which  generate  the  surfaces  0  =   const. 
3o   The  vectors   a*^   satisfy  the  linear  algebraic 
equations  (i|o20)o 
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5 .   Generalized  P.ropi;resslng  Wave  Solutions  of  the  Cauchy 
Problem  for  Systems  of  Higher  Order  with  Multiple 
Characteristic  So 

In  this  section  we  shall  use  the  formal  construction 
of  Section  l\.   to  obtain  the  solution  of  a  Cauchy  problem 
Involving  an  arbitrary  functlono   We  shall  see  how  this 
solution  can  be  carried  out  for  a  class  of  hyperbolic 
equations  characterized  by  an  algebraic  condition. 

¥e  label  these  equations  "equations  of  type  H." 
These  are  equations  of  the  form 

(5ol)         ^u  E  A(D}  u  +  B(D)  u  =  0  5 

where  the  principal  part  A(D)   satisfies? 

lo  The  hyperplanes  x     ~   const «   are  '"spacelike  o " 
We  interpret  this  condition  as  meaning  that  the 
determinant  of  coefficients  In  equation  (Soil) 
does  not  vanish. 

2.  The  system  (I|.«13)  reduces  to  ordinary  differen- 
tial equatlonso— / 

3o  The  multiplicity  of  each  root  of  the  characteristic 
equation  is  constant. 

I4..  All  coefficients  are   C°°   functions  of  x. 

The  following  types  of  equations  are  included  in  class 

Hs 

97 — -— — 

-'  This  condition  is  satisfied  if  A    is  proportional  to 
the  identity.  ^o 
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(a)   A   first-order   system,    possibly  with  multiple 
characteristics,   where   for   fixed      (^n , « « •  ,i^    ),      the    set 


of  nullvectors   R 


a(3 


(p  =  l,oa,,r  ,  a  =  1,.,.,M)   of 


the  characteristic  matrices   A(D0  )   form  a  complete  set. 

(b)  A  single  equation  of  order  m  with  m  real  and 
distinct  characteristics  for  fixed   (^-joo.,]?^  ) . 

As  before,  we  denote  x   by  t   and   (x  ,...,x  ) 
by  y.  Consider  the  Cauchy  data 

00 

f.((2((v))  t 


iB.2) 


u(o,y)  =  :2_  f.(j2^(y))  b^  , 
j=0  J 


00 


u^(0,y)  =  ;gl  fj_^  (^(y))  b^  , 


00 


V 


u^„.l(0,y)   =  S  ^,J-m*l'^<y"  "m-l 


where  the  following  additional  conditions  are  satisfied: 
So   f  (s)   is  an  arbitrary  function,  and 


&^J^^) 


^j_l(s) 


( J    =   -m+1 , . , o ) . 


,00 


60      0      is   a     C"^       function   of     y,      and    the   gradient 
of     0     with  respect   to      y     does   not   vanish. 

7.      The   vectors     bP(y)      are      C°°      functions    of     y. 

As   in  Section  3,    we   first    construct   all  possible 
solutions   of   the    nonlinear  first-order  partial    differential 
equations 
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(5.3) 


det  {A{J)0^))      =   0 


(a  =  1, . . . ,M), 


with  the  initial  condition 
(5.i+)         ^''(0,y)   =  0ij) 


(a  =  1, . . .,M) 


This  is  equivalent  to  solving  the  associated  equations 
homogeneous  of  degree  1  (see  equation  (I|.»15))! 

(5.3')         0^  +  X''{0^,,.,J^)      =   0      (a=l,...,M), 

or   the    associated   canonical    system   of   ordinary   differ- 
ential  equations   along  the   rays: 


dt     ^     ,        ^     ^     ^ 
dr    -        ^      dr      -     ^^a 


(5.5)^ 


dp° 


(li  =  1, ..  .,n) 


( V   =   0,1, .. . ,n) 


(a=l,  ...,M), 


with  the  initial  conditions 

f   t(0)   =   0  ,   x^'(O)   =  y^" 


^ 


(5.6) 


1 


f  ( tx— 1 ,  • .  o ,  n ) . 


¥e   assT;mie   that  there    is   a  number      t,    >   0      such  that   the 
solutions   of    (5.3)    are   regular   for      0  <   t   <  t, .      Such  a 
number    always   exists    if   the    Gauchy   data    {S»2)    has    compact 
support. 

Using   the   techniques   of  Section  ii,    we    can  construct 
a  formal    solution   of   the    Gauchy  problem.      We    assijme   that 
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w(x)      has   the   form: 
(5.7)  w(x)      = 


YZ  f.(^    (x))    a'^J    . 
j=0      a=l      J 


To  each  of    the    functions      0    (x)      are   associated  vectors 


,ap 


and     L 


ay 


such  that 


A(D^^)    R°-^     =0        (p  =  l,...,r^;    a  =  1,...,M), 
L'^^  A(Dj2(°)    =   0  (y  =   l,.o,,r^;    a  =   1,...,M). 


We   require   that     a      (x)      have    the   f 

r 
(5.9) 

Then 


orm 


ao 


(x)      =     JZ  ^R^'^^x)   R'''^(x)      h     o^^°  R 


ap 


p=l 


w(0,7)   =  f    (^(y))[7~cr„^°R^P]-t-T~  f,(^)[...]. 


^P 


J=l 


M  00 

wiO,y)    =  f   .  (^(y))[y-(-A^)gr;^°R^P]+J-  f .  ,(0)[...   J, 
^  "-^  a=l  P  .1=1     J"-^ 


(5.10). 


M 
^  m-l(°»y^    =  ^W^(7))    [Z:   (-^^)"^-Vp'^°R^P]    + 


00 


V 


-^  z:  i'j+i-^(i2')  [...J . 


Identifying   the   first   coefficients    in  equations    {B»2)    and 
(5olO),    we    obtain  the    system   of   equations: 
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M 

T- ClO/^       \    T3C13      _      i_o 

X_  cTp      (0,7)   R   ^     =     b 

a=l      P  ° 


■^r —  ao/    -.a\    „as     _     ,0 

a=l      ^  -^ 


(sai) 


M 
J-  o^^^°(-X^)™-l   R^P      =     b°   ,     . 

This  is  a  system  of  m  vector  equations  with  k 
components  each  for  the  m   k   scalars   cr;   (0,y)<,   By 
hypothesis,  the  determinant  of  the  system  does  not  vanish. 
Now  we  shall  verify  that  this  hypothesis  is  satisfied  in 
cases  (a)  and  (b). 

(a)  In  case  (a),   m  =  1.   System  (5«11)  consists  of 

a  single  vector  equation.   There  are  just   k  vectors  R  ^, 
and  by  hypothesis  they  form  a  complete  set.   Hence  the 
equations  have  a  unique  solution. 

(b)  In  case  (b),   k  =  1„   System  (5oll)  consists  of 
m  equations  for  the  numbers  cr  ^^ o      The  determinant  of 
the  system  is  the  Van  der  Monde  determinant  of  the  numbers 
"k   joooyX  „   By  hypothesis,  these  numbers  are  distinct,  and 
hence  the  determinant  does  not  vanish. 

Hence,  in  each  of  these  cases,  the  scalars   or    can 
be  determined  from  (Soil). 

Now  the  factors  (T„  can  be  determined  along  the 

rays  which  start  at   (Oj,y)   by  integrating  the  systems  of 
equations s 
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dcr: 
dr 


ao 


(i,.l8)   p^^P  :ZL.p^TP^^ao  ^  ,   (Y  =  l,...,r^) 

(a  =  1, . .  ,,M). 


Then  (5.9)  yields   a'^^Cx).   This  choice  of  a^°   insures 
the  compatibility  of  the  systems; 


([1..11J  AlDjZl^ja^^+A^CD^^joa^^+Cl^^ja^^  =  0,  (a  =  1,...,M), 

P       -^ 

and   a    can  be  determined  modulo  a  right  nullvector  of 
A(D^^): 


al     _     ,  al    ,    _  al  ^a-P 
L  -     h       ■*"   ^R       R        J 


(A  =  1, .. .,M). 


al 


Now  we    can  obtain  the   initial   conditions   for      o_  by 

identifying  the    coefficients   of   the    second  terras   in    {5<>2) 
and   (5«10); 


a=l      ^ 


M 


IZ  o-"^  R°P  =  bi  -  2Z  h"^  , 


a=l 


\     M  ^             TO               T  M        X    ao                         , 

ZI  (->^    )    a^*"^   R''^      =      b,l  -  IZ  (^  +    (-A'^)h^), 

J  a=l  ^                           ^  a=l       <^^ 

(5.12)/ 


t=l  P  "^-^        ^  ^^ 


V 


+    (-A^)^-l   h^^J    . 


ttj 


In  general,    the    Initial    conditions   for      en    ''      are 
determined  by   solving  a   system   of  equations   of   the    form: 
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M  M 

M  ■■  M 

a=l  ^^  -^        a=l        -^ 


(5.13V' 


where      S    ^      depends   linearly  on     a      ,..,j,a     ^  and   their 

derivatives  with  respect   to     to 

The   factors      en    ^      can  be    determined  by  integrating 

the    systems   of   equations? 

d  "■  °'^ 
(i^,21)      p^yP     _1l_  +  p^YP   ^aj   ^   _^Y(A    ,h^   +  chJ   +  K^j, 


(y^Ijo../"    ),    (a=l,.o.,M). 


Then  the  equations 


{ho20)  A(D^'')  a''-5"^^+  A  «>  Da""^'  +  Ca^'J  +  k'^J   =  0  , 

can  be  solved  for   h  -^     ,      Thus  the  terms  of  the  series 
(5o6)  can  be  constructed  in  succession. 

The  argument  of  this  section  has  shown  that  if  we 
assign  Cauchy  data  in  the  form  of  a  wave  to  a  hyperbolic 
system  of  type  H,  and  if  the  formal  solution  is  assiimed  to 
consist  of  a  sum  of  generalized  progressing  waves,  then  the 
formal  solution  has  a  unique  decomposition  into  a  sum  of 
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such  waves.   Each  of  these  waves  has  a  family  of  charac- 
teristic surfaces  of  the  system  for  phase  surfaces.   We 
shall  apply  this  formal  construction  in  the  next  section. 
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6,   Applications  of  the  Formal  Construction 

A.  Propiressinp:  waves  and  propyesslng  waves  of  higher  degree. 

It  can  happen  that,  by  a  fortunate  combination  of 
circixmstances,  all  but  a  finite  number  of  the  distortion 
coefficients   a^'*^   vanisho—''   In  this  case  the  formal 
solution  is  an  actual  solution,  or  rather  a  class  of 
actual  solutions,  since  the  wave  form  f^^s)  ^^J   tie 
chosen  arbitrarily.   Such  a  class  of  solutions  is  called 
a  progressing  wave  of  higher  degree*   If  only  the  first 
terms   a^'^   are  different  from  zero,  the  solutions  are 
called  progressing  waves.   Examples  of  progressing  waves 
and  of  progressing  waves  of  higher  degree  were  given  in 
the  introduction. 

B.  Power  series  expansions. 

1   1 
If  we  choose  f  ,(s)  =  -r-r  s-Jj,   then  our  formal  expansion 

becomes  a  power  series  in  the  functions  0  (x).   In  many 

cases  this  series  converges  (see  Section  7  )<!  and  thus 

represents  an  actual  solution  of  the  Cauchy  problem.   In 

this  case  our  method  of  solution  is  similar  to  that  used 

by  J.  Hadamard  [5 J  to  construct  the  fundamental  solution 

for  equations  of  second  ordePo 

C.  Asymptotic  expansions. 

If  we    choose      fAs)   =    (l  .O'^^e         ,      then   o\xr    formal 

— /   Such  a   situation  has  been  analyzed  for   the  wave 
equation  by  P.    Go   Friedlander    [6], 


expansion  becomes  an  asymptotic  expansion  of  the  solution 
of  the  Cauchy  problem  in  terms  of  the  frequency  parameter 
^-^o   In  this  case  our  method  of  solution  of  the  Cauchy 
problem  is  identical  with  that  used  by  P«  D.  Lax  [10]  to 
obtain  such  asymptotic  expansions  for  solutions  of  systems 
of  first-order  equations, 

D,  Weak  solutions  of  problems  with  distribution  initial  data. 

Since  the  distortion  factors   a  ^^   do  not  depend  upon 

the  wave  form  f  (s)   involved  in  the  initial  data,  we  can 

o  ' 

allow  generalized  functions   f  j,   which  are  limits  of  smooth 
functions  in  a  topology  which  respects  addition  and  multi- 
plication by  a  smooth  function.   In  particular,  if  f    is 
a  distribution,  then  the  corresponding  formal  solution  of 
the  Cauchy  problem  enables  us  to  construct  a  weak  solution 

of  the  Cauchy  problem. 

P   P 
More  specifically,  let  f  (s)  =  d  /ds   g(s)   where 

g(s)   is  merely  a  continuous  function,  and  P  is  a  positive 

integer.   We  assume  further  that  only  a  finite  number  of 

the  factors  bv   are  different  from  zero  (see  eq,  (5«2))o 

Let 

j=0   a=l      J 
The    factors     a   ^      have  been  constructed   in  such  a  way   that 

^'(u-^)    =  2Z  ff  T  n,+l  (^'')[A^»Da^'^+Ca'''^+K^'^j+  J-     f  .  (^2^^)  [,„  ]]. 
^  [^   J-m+1  p  ^=jri+2    -  J 

M     /~  J-P  "1 

=  .E:fgj.p„^-Hi(i2^'')tA  "Da'^'^+Ca'^'^-HK'^J-H  ^Z       g/^'')[...]  \. 
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Thus     (u  )   has   J-P-m+l  derivatives  in  the  ordinary 
sense.   We  now  define  a  function  v(x)   by  the  conditions 

(6„1)  /(v)   =  "   .^(u*^)  , 

(6.2)      v(0,y)   =  v^(0,y)  =  o..  v  ^^i^iOsJ^      =  0  • 

This  is  a  Cauchy  problem  and  it  has  smooth  data  provided 
that   J   is  taken  sufficiently  large.   Therefore  the 
standard  existence  and  uniqueness  theory  (see  Leray  [11  J) 
implies  that   v  exists,  is  uniquely  determined,  and  has 

as  many  derivatives  as  we  like  if   J   is  chosen  suffici- 

11/ 
ently  large o — ^   Now  we  define   u(x)   by 

u(x)   =  u  (x)  +  v(x)  , 

This  function  u(x)   is  the  unique  solution  of  the  Cauchy 
problem  (^d},  {5»2)o 

We  can  summarize  our  results  in  the  following 
Theorems   If  the  singular  part  of  the  Cauchy  data  has  the 
form  of  a  wave  (So 2.)^    and  if  there  are  no  caustics  (which 
is  the  case  for   t  <  t,  ),   then 

lo  The  singular  part  of  the  solution  has  the  form  of 
a  sum  of  generalized  progressing  waves.  If  j?(y)  is  the 
phase  function  for  the  initial  data,  then  the  phase 

— ^    Once   J  has  been  chosen,  the  differentiability 

conditions  on  the  coefficients  of  the  equations  and 
on  the  initial  data  can  be  relaxed.   We  shall  not 
attempt  to  show  how  many  aerivatives  are  actually 
required. 
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functions  for  the  singular  part  of  the  solution  are 

solutions  of  the  characteristic  equation  of  the  system 

which  agree  with  0{j)      for  t  =  Oo   Thus  if  the  family 

of  curves   C(c)   is  defined  by  the  equation  0{j)   =   c, 

then  the  solution  is  a  smooth  function  except  along  the 

characteristic  surfaces   S  (c)  which  emanate  from  C(c). 

2o   The  distribution  f  (s)   which  gives  the  wave 

o  ^ 

form  of  the  data  along  C(c)   also  gives  the  wave  form 
of  the  solution  along  S  (c}o   Consequently 

3o   The  solution  has  the  same  smoothness  properties 
as  the  initial  data. 

[|.o   The  distortion  coefficients   a  ^   are  determined 
by  solving  linear  first-order  ordinary  differential 
equations  along  the  rays  which  generate  the  characteristic 
surfaces   S  (c). 

Moreover, 

^o   These  coefficients  are  independent  of  the 
particular  wave  form  in  the  Initial  data. 

6„   The  coefficients   a*^*^   are  identical  with 
coefficients  which  appear  in  the  asymptotic  expansion  of 
the  solution  of  a  Gauchy  problem  with  rapidly  oscillating 
Initial  datao 
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7o   Proof  of  Convergence  In  the  Analytic  Case. 

¥e  have  pointed  out  In  the  introduction  that  pro- 
gressing vrave  expansions  are  meaningful  regardless  of 
whether  or  not  they  are  convergent.   Nevertheless,  the 
situation  is  simpler  when  such  expansions  converge.   It 
is  apparent  that  convergence  of  such  expansions  must  be 
connected  with  analytlclty  of  the  equation  and  the  data, 
since  a  power  series  expansion  can  be  obtained  from  a 
generalized  progressing  wave  by  appropriate  choice  of 
the  wave  form.   In  fact,  we  can  show  that  if  the  coeffi- 
cients of  the  equation  are  analytic,  and  if  the  distortion 
coefficients  of  the  initial  data  and  the  initial  values  of 
the  phase  function  are  analytic,  and  if  the  power  series 
corresponding  to  the  initial  values  of  the  wave  is  conver- 
gent, then  the  wave  converges  in  an  appropriate  neighbor- 
hood of  the  Initial  surface.   Note  that  no  assumptions 
need  be  made  about  the  wave  form.   In  particular,  the 
theorem  applies  to  progressing  wave  expansions  where  the 
wave  form  is  a  distribution.   The  theorem  also  applies  to 
equations  which  are  not  hyperbolic,  except  that  for  non- 
hyperbolic  equations  the  Cauchy  data  cannot  be  chosen 
arbitrarily.   The  proof  consists  in  transforming  our 
problem  into  an  analytic  characteristic  initial  value 
problem,  and  applying  an  extension  of  the  Cauchy- 
Kowalewski  theorem.   For  simplicity,  we  present  the  proof 
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for  a  first-order  system. 

First  we  shall  give  an  extension  of  the  Cauchy- 

Kowalewskl  theorem  to  apply  to  the  characteristic  initial 

value  problem.-—/   In  order  to  clarify  the  presentation, 

we  shall  first  prove  a  theorem  which  is  weaker  than  the 

one  we  need. 

Theorem  7.1.   Let   A^(x)   (v  =  0,..,,n),   B(x)   and  C(x) 

be   k  X  k  matrices  whose  elements  are  analytic  functions 

of   (x°,...pX^)   in  a  domain  D  in   (z,x) -space  which 

includes  the  origin.   Let   C(x)   have  rank  k-1,   and  let 

R(x)   and  L(x)   be  unit  right  and  left  nullvectors  of 

C(x)   respectively,   ¥e  assume  that  the  scalar  LA  R  7^  Oo 

We  may  also  assume  without  loss  of  generality  that 

n  -\ 

y^  La"^R  >  0  in  D.   Let  b(z,x  ,..,,x^)   be  a  vector 
v=0 
with  k   components,  each  of  which  is  analytic  in  D, 

with  b(0,x  ,...,x  )  =  Oo 

Consider   the    system   of      k     linear   first-order  partial 

differential  equations 

(7.1)  c(x)   ^4l^     =     Z:  A^(x)  i^ii^^  +  B(x)   u(z,x), 

'^^  v=0  ^x 

with  the  subsidiary  conditions 

(7.2)  u(0,x)   =  0, 

(7.3)  [u(z,0,x^,,..,x'^)  -  b(z,x^,.o.,x'')J  •  R(x)   =  0. 


— /   As   far   as   I   know,    such  a  theorem  has  been  proved   only 
for    second-order  equations    (see  Hadamard    [7 J   PP.    77-81) 
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If  the  above  hypotheses  are  satisfied,  then  (7»1),  (7«2) 
and  (7»3)  have  a  unique  analytic  solution  in  an  appropriate 
neighborhood  of  the  origin. 

Proof:   First  we  remark  that  (7.1)  and  (7»2)  consti- 
tute a  characteristic  Cauchy  problem.   Now  we  shall  show 
that  equation  (7«3)  removes  the  inde termination  in  the 
formal  expansion  of  the  solution. 

It  follows  from  (7.1)  that 

il.k)  L  (  y^  A^  ^  +  Bu)  =  0  . 

v=0    jx 

If  (7«4)  is  satisfied,  then  ^n/r)z      is  determined  in 

terms  of   u  and  c)u/^x   modulo  a  right  nullvector  of 

C.   We  write 

(7.5)  u(z,x)   =  h(z,x)  +  o-(z,x)  R(x)  , 

where   o-  is  a  scalar  factor.   Prom  (7.1),  we  may  obtain 

(7o6)  h^  =  o^^(h)  +  ,^2^°"^  ' 

where     'X_         and     ,<_  „      are   first-order   partial   differential 
operators  which   Involve    only  the  x-variables.      Equation 
{7ol\.)   may  be  rewritten  in  the   form 

(7.7)  ff  +  P(x)   o-    =    c/3(h)    , 


?       LA^R,,    ,      and     c^, 
is   a  first-order  partial   differential   operator  which 


where      dxVd'^     =  LA^R,      P(x)   =  ^^  ^^^    >      and     c^, 

v=0  ^ 
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involves   only  the  x-variables.      Eqixatlon   (7.3)  may  be 
rewritten  in  the   form 

(7.8)  cr  (z,0,x-'-,...,x^)     =      (b-h)  •  R   . 

The  formal  expansion  of  u   can  be  computed  from 

(7.6),  (7.7)  and  (7.8).   It  remains  to  show  that  this 

formal  expansion  converges  in  a  neighborhood  of  the 

origin.   To  do  this,  we  shall  construct  majorants  for 

h  and  o" . 

J3- 


Let   y  =  >   X   o   Let  M-,/(l-y/p,  )   majorize  the 
v=0  V         / 

norms  of  the  coefficients  of  ~^ ,   and  ^p.   If  S 

majorizes  <y ,      and  if   H   satisfies 

(7.10)  H(0,x)   =  0  , 

then   |hl  <  H.   In  other  words,  if  we  have  a  majorant 
for   o  ,   then  we  can  easily  obtain  a  majorant  for   h. 
Note  that  in  order  to  estimate  the  coefficients  of   z 
in  the  expansion  for   h,   it  is  only  necessary  to 
estimate  the  coefficients  of   z  ,  o  o . ,  z '"    in  the 
expansion  of   cr  o   Note  also  that  if   S   depends  only 
upon  z   and  y,   then  H  depenas  only  upon  z   and  y. 

Now  we  proceed  to  estimate   o-  in  temis  of   h. 
From  (7.7)  and  (7.8),  we  may  represent   o   in  the  form 
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(7.11)      Or-(z,x)   = 


Jk(x,t-}:^[h(z,x(Tr))]  dr  + 


+  m(x)  [b(z,x' )  -  h(z,x' ) ]  . 


Here   k(x,c)   and  m(x)   are  bounded  functions,  and  x' 

is  the  intersection  of  the  ray  containing  x  with  the 

surface  x  =  Oo   Since  the  vectors   (LA  R,...,LA  R)   and 

(!,.«,, 1)   are  not  orthogonal  (see  the  hypothesis)  we  may 

introduce   y"  =  ^   x  ( r)   as  variable  of  integration 

v=0 
instead  of  f  .   It  follows  that  if   H(z,x)   is  a  majorant 

for  h,   H  is  a  function  of   z   and  y   only,  and  M. ,  p. 

are  properly  chosen,  then 


(7.12)    |o-  I    < 


M, 


n 


7^  (^  ^^ 


1-  nil    ^j=o 
P2 


1.  ml 


H 


dy"  + 


M 


k 


1-  ^ 


H(z,x')    + 


zM 


5 


(1-  -f)(l-  -f) 
P5  P6 


But 


>   H  =  (n+1)  -T—  Ho   Hence,  integrating  by  parts. 


v=0 


(7.13)  kl  <  -  J 


(n+l)M,  1^(1-  -zbr)  - 


'2  ^^'   ^ 
P2 


M 


1. 


P3 


H(z,x(y"))dy"+ 


M2(n+1) 


M2(n+1) 


M 


k 


+    p^  H(z,x)=  \^,i^\   H(z,x')+  — ^  H(z,x')  + 
P2  P2  PI4. 


zM 


£ 


(1-  ~)(i-  -f ) 

P5     P6 
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If  Mp   is  chosen  large  enough,  and  pp   is  chosen  small 
enough,  then  the  integrand  in  (7«13)  can  be  made  positive, 
and  the  term  involving  H(z,jc' )   can  be  made  negative. 
Hence  the  inequality  is  improved  if  these  terms  are 
dropped.   Thus  we  obtain  an  inequality  of  the  form 


(7.II4.) 


M, 


icr 


1-  ^ 


H(z,x)  + 


zM, 


(1-  :^)(1-  ■^) 


We  may  set 
(7.15)    S  = 


M 


1-  ^ 


^  H(z,x)  + 


zM 


L 


(1-  -^}(i-  ■^) 


'7 


'7 


,^ 


Note  that  the  coefficient  of  7.^      in  the  expansion  of   cr 
may  be  estimated  in  terms  of  the  coefficient  of   z' '  in 
the  expansion  of   H,   Also,  if   H  depends  only  on  z 
and  y,   then   S   depends  only  on  z   and  y. 

We  assert  that  a  majorant  of   h  is  given  by  an 
equation  of  the  form 


(7.16) 

with 
(7ol7) 


H. 


M 


1-  ^ 
P 


n 


.)H 


+  H  + 


V=0  r)x 


1-  ^ 
PJ 


H(0,x)   =   0 


A  majorant  for  o-  then  is  given  by  (7.15).  The  constants 
M  and  p  are  obtained  by  a  formal  substitution  of  (7.15) 
into    (7.9).      The   proof   of  this    assertion  proceeds  by  means 
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/, 


of  an  induction  on 

Assume  that  the  expansions 


jjv.c  -^   =  ^  Hj(x)  z 


and 


majorize  the  first  -^-1   terms  in  the  expansion  of   h 


and  o^o   Ass'ume  further  that  H^'*"    was  obtained  from 
(7.16)  and  S^"^"-"-^  was  then  obtained  from  (7.15).   We 

will  show  that  the  corresponding  statements  are  true  of 

(O        (^-) 
H     and  S   .   It  follows  from  the  remark  following 

(7.9)    that     H^'^^     may  be   obtained  by  substituting   (7.15) 

iJ) 
into  (7.9),  i.e.   H^^''   is  given  by  (7.16).   In  a  similar 

way,  we  see  that  S^   is  given  by  (7.15).   It  only 

remains  to  show  that  the  induction  can  be  started.   This 

is  indeed  the  case,  since   H    =  0  and  S    =  M_z  . 

Since  (7.16),  (7.17)  obviously  have  a  unique  solu- 
tion, we  conclude  that  the  expansions  of  h  and  cr 
converge  in  a  sviff iciently  small  neighborhood  of  the 
origin. 

As  we  have  mentioned.  Theorem  7*1  is  not  strong 
enough  for  our  purposes:   we  shall  have  to  deal  with  the 
case   where   the  initial  manifold  is  multiply  character- 
istic.  In  order  to  avoid  clumsy  algebraic  hypotheses,  we 
shall  assume  that  our  equations  are  already  in  a  form 
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corresponding  to  (7»6),  (7.7)  and  (7.8)0 
The orem  7»2s   Consider  the  set  of  equations 


1       2 


(7.18)  h^   =   -  j_(h,.  +   :  2((^  ,.o.,cr^)  , 

{7ol9)    ^^^  +  F^'U)    o-°-     =  ,.-.  ?(h)       (a=l,...,r), 

(7.20)    o-°(z,0,x-^,.,.,x^)   =  E'^(x)h  +  F'^(x)b  ,  (a=l,..,»r) 


Here     h     is   a  vector,    and     c         may  be   either   a  vector  or 

a   scalar.  -,,      -p      and  ^      are   first-order  differential 

operators  which  involve    only  the   x-variables,      P°'(x),    E°'(x) 

and     F    (x)      are    coefficients  which  may  be    scalars,    vectors 

or  matrices,    depending  upon  whether     o  is    a  scalar   or   a 

vector.      All   functions   are   assumed   to   be    analytic   in  a 

neighborhood   of   the    origin  of      (z,x) -space.      We    assume 

,    o 
further   that     — ■, — r  ^  0.      By  means   of   an  affine    transfor- 

^'^  °'  -IL.  V 

mation,  we  may  insure  that   y  =  2. ^   ™^7  ^®  introduced 

v=0  /  ^ 

as  a  curve  parameter  in  (7.19)  instead  of   ~    .   If  all 

of  the  above  hypotheses  are  satisfied,  then  equations 
(7.18),  (7.19),  and  (7.20)  (and  of  course  the  correspond- 
ing Cauchy  problem)  have  a  unique  analytic  solution  . 

The  proof  is  an  obvious  modification  of  the  proof  of 
Theorem  7.1,  and  will  be  omitted. 

After  these  preliminaries,  we  can  easily  prove 
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Theorem  7o3"   If   Xu  =0   is  a  system  of  linear  first- 
order  hyperbolic  equations  with  analytic  coefficients 
of  type   H   (see  Section  5),  and  Gauchy  data  is  prescribed 

"1  Tn 

on  the  spacelike  manifold   t  =  0   in   (t,x  , .oo,x  )   space 
in  the  form 


CO 


(7o21j         u{0,x)   =  X^bhx)    f  Cpixj)  , 

j=0        ^ 


where 

le   f    is  the  m-th  derivative  of  a  continuous  function, 
o 

3»   f  .(0)  =0   if   j  >  m, 

ko      p(x)      is    analytic  J    and  grad      0^0    , 

i  Tr22-      i      z^ 

So      b'^lxj      are    analytic,    and   the    series      >       b'^ixJ-r-r- 

converges, 

then  the    solution  has   the   form 

k        00  . 

(7022)  u(t,x)      =     T~  ^~   a^J(t,x)    f,(0^(t.xj) 

a=l   j=0  J 

where   j2  (tjX.)   are  the   k   solutions  of  the  characteristic 
equation,  with  0  (0,x)  =  i3(x),   and   a  ^(t,x)   are  coeffici- 
ents which  are  independent  of   f  (0)  o 

Proof;  We  first  remark  that  it  is  sufficient  to  prove 
the  convergence  of  the  series 

(7023)  v(t,x)   =  y^  T~   i^-  a^J(t,x)  , 

a=l  j=0    J° 

since 
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ilcZk)  \f^(0}\      <  M  M__j_  ,        for   j  >  m  . 

a^^      are  determined  from  the  following  set  of  equations: 


(7.2M         A^^^  a^J^l   =  -^{B.^'h  (j  =  -1,...), 

Where     A^"^''    =  A^^°'   ,         a'^'^     must   also   satisfy  the    initial 

condition 

k 
(7o26)  ZZ  a°'J"(0,x)      =     b^  (j   =  0,l,oo.). 

a=l 

In  order   to    deal  with  this   infinite   set   of  equations, 

we   introduce    the   generating  fionctions     g    (z|tj,x)      and 

b(z,x) s 

_ro_  j+1 

(7.27)  g^{z^t,x)      =     2„  a'^Jlt.x)  77+5-71  » 

(7<.28)  b(z|x)      =     ^>_  bJ(x)  -rrmT   " 

The    infinite    system    (7-25)    is    equivalent   to    the    set 
of     k     partial   differential   equations s 

(7.29)  A^^'^  4f^     =      -c£ig°-)  (a  =  l,...,k), 

and  the  initial  conditions  become 

k 
(7o30)  ZZ   g  (zpO,x)  =  b(z,x)  . 

a=l 

By  definition,   g'^  satisfies  the  condition 
(7.31)  g^(Oi;t,x)   =  0  (a  =  l,..o,k) 
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For  fixed   a,   equations  (7«29)  and  (7.31)  constitute 
a  Cauchy  problem  for   gCzjtjX),   where  the  initial  mani- 
fold Is  given  by   z  =  0,   not   t  =  0.   Since  0   satisfies 
the  char'acteristic  equation,  the  matrix  A     has  a 
vanishing  determinant,  and  hence  the  manifold  z  =  0   is 
characteristiCo   Hence   g  (z;t,x)   need  not  vanish 
identically,  even  though  it  must  satisfy  a  homogeneous 
equation  and  assume  zero  initial  values  (for   z  =  O). 
Indeed,   g   is  not  uniquely  determined  unless  conditions 
are  prescribed  on  another  surface,  in  this  case  the  mani- 
fold t  =  0„   This  additional  condition  (7.30)  links  the 
systems  (7«29). 

In  order  to  treat  condition  (7«30),  it  is  expedient 
to  write  our  equations  in  still  another  fonrio   Let  G  be 


a  vector  with 


components: 


(A 


2     2 
let   C  be  the   k  x  k   matrix 


C  = 


¥e  have 


/  a'i' 


0 
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o4i  = 

ciz 

"    ^(G)    , 

GlOjtjX) 

=      0      , 

a=l 

=     b(z|x) 

(7.32) 

(7. 33) 
(7»3ij-) 


This  is  another  characteristic  Inital  value  problem, 
but  in  this  case  the  manifold  z  =  0   is  characteristic 
with  multiplicity  k„   Equation  (7»3l4-)  gives  just  k 
additional  conditions  on  the  manifold   t  =  0„   The 
hypotheses  of  the  theorem  imply  that  the  coefficients 
in  equation  (7o29)  and  the  data  (7»3^)  are  analytic* 

Since  our  equations  are  of  type  (H),   the  coefficients 
a'^'^   are  determined  from  equations  ([4.0 20)  and  ( 1^.021  )o 
These  equations  correspond  exactly  to  (7ol8)  and  (7.19), 
while  our  condition  (7o30)  immediately  yields  an 
equation  corresponding  to  <7o20)o   Thus  we  can  apply 
Theorem  7o2,  and  the  series  (7«23)  for   v(t,x)   converges 
if  0'^fo  o ,  ^0^^t,    and  x  are  sufficiently  small. 
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8.   The  Rlemarm  Function. 

In  this  section  we  shall  give  a  method  for  construct- 
ing the  Riemann  function  for  equations  of  type   H. 

We  may  summarize  our  results  as  follows  J 

1.   The  singular  part  of  the  Riemann  function  can  be 
represented  as  a  superposition  of  waves  whose  phase 
siirfaces  are  regular  characteristic  surfaces. 

2.,   In  the  analytic  case,  the  Riemann  function  can 
be  represented  as  a  superposition  of  convergent  progress- 
ing waves.   These  convergent  progressing  waves  are 
constructed  in  a  manner  similar  to  Hadaraard's  method  for 
equations  of  second  ordero 

3.   The  Riemann  function  is  a  smooth  function  except 
if  its  argximents  lie  on  the  same   ray. 

A.   Differential  operators  and  their  adjoints. 

We  shall  make  use  of  the  notation  of  Section  l\..      Let 

(8.1)  Q(p)   =  JZ  ^J^^   P""  . 

a 

where      q^(x)      are      k  x   k     matrices,       q      =0      if       lal    >  m 
and     a  =1.      Let 

TTl  ,  O  5  O  ,  o  »  .  ,  O 

(8.2)  u      =     Q(D}u      =     2Z  q    (x)    D%   . 


We   define       ;^     v      by 


a 


-r  •«>        _     ■'T —    /    T^  \  a 


(8o3)  .^    V     =     JZ   (-D)      V  q^ 

a 
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.^^      is  called  the  formal  adjoint  to  -yC  «   >.'_  and  <__ 
satisfy  the  relation 

(8.1;)   vc^u  -  (w*v)u  =     Yl^'    P^(v,u)  , 

i=0  ^ 

where      P    (v,u)      is   a  bilinear   differential   operator 

whose    coefficients    depend   on     q   (x). 

The  right-hand   side   of    (Soi^.)    is    clearly   a   divergence. 

Let      V     be    a  region  of  x-space  with  boundary      So      Gauss' 

theorem  implies   that 

(8.5)  /     [v/u-(  :<£%)u]    dx     =         (     5~"  P^(v,u)    dS,    . 
-V  S   i=0  ^ 

This   is    the   fundamental  relation  between  a   differential 

operator    and   its   adjoint,    sometimes    called   Green's 

identity.      If      yikiy ^,^)    =  v^(  z^'^^y^   =   ...    =  ^  ^,2(^2*^^    ^  °' 
and     V      has    compact    support   in  every  plane      t   =  const., 
then,    integrating    (8.I4.)    over   the   region         ",    <      t     <     C'^ 
we  have 

(8.6)  /  [v-(u-(-^*v)u]dx  =      \{-\?'-\  ^   n(~P,7)u('rp,y)dy 
'^llt<C^2  "^ 

-  yP°[v(',^^,y),u(-r^,y)J   dy, 

since      a  =1. 

Tn,  o,  „  o  .  p  o 

B.   The  Riemann  function  and  its  construction. 

To  any  hyperbolic  system,  we  can  associate  a  funda- 
mental solution  which  bears  a  particularly  useful  relationship 
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to  the  Cauchy  problem,  namely  the  Riemarm  functlorio   The 
Riemann  function  corresponding  to  the  operator  -^.^ 
satisfies  the  following  conditions: 


■5», 


if  t  <  r'. 


(8„7) 


k(t,y,ei>j/)  =  0 
k(  ,y,^i>/)  =  k^(^;y,'',^i)  =  .00  =  ^^^-z^"^''^'^'^]^"^' 


ra-1. 


3'  \m-l^^>^'^''^   =  (-l)^'^"^6(y-.)l 


ko      k(t,y,;,n)  =  0 


if   t  >  ^  . 


5(y}   denotes,  of  course,  the  n-dimensional  Dirac  distri- 
butiono   The  uniqueness  theory  for  hyperbolic  equations 
implies  that   k  has  compact  support  in  each  plane 
t  =  consto   Inserting  (8o7)  into  (8o6)„  we  have 


(8„8)   u(  r^,'/)  =   /' 


"  a<t^^ 


k(t,y. 


)  'u  dt  dy  + 


+  yp°[k(c-^,y,r2,>/)  *  ^^"i^y^J  ^y 


It  is  because  of  this  formula,  which  expresses  any 
function  u(x)   in  terms  of   '.u  and  the  Cauchy  data 
of   u  on  the  hyperplane   t  =  ?%  ,   that  the  Riemann 
function  is  so  importanto 

We  may  define  the  Riemann  function  for  the  operator 
oC-   as  follows s 
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r 


(8 


.9)S 


i.Xk^(t,7,-,^)  =  0  if  t  >  r', 

2,  k^*(cy,r,)^)  =  k*(r,y,r,>|)  =  ..»  =  k*  ^r,y,'  ,)|)=0, 

3.  k*j^_^(r,y,(  ,/•/)  =  6(y-ri)  I  , 
i;o   k(t,y/;,u)  =0  if  t  <  l\ 


We  set  v(t,y)  =  k(  t,y,-C2»'',2  ^  '   u(t,y)  =  k  (t,y,  r]_,>^) , 
and  substitute  into  (8o6)c   This  yields  immediately  that 


(BolO) 


k'*(r-2,^i2'''i*':'l^  "  ^^'^l*''/l-'^2*'?2^ 


This    symmetry  relation  exhibits   the    close   relationship 
between   a   system  of   equations   and   its   adjoint^    and 
demonstrates   that   the    function     k(t,y,~, M)      satisfies 


(Boll) 


^^^^k(t,x,::,!^) 


0 


(^  >      t) 


with  respect   to   the   variables      ■>       and     >■/    » 

The   Riemann  function  can  be    constructed  by  means    of 
progressing  waves,    provided  that   the    5-function  can  be 

represented   as    a   superposition   of  waves  with  regular  phase 

13/ 
surfaces,    for  example   plane  waveso-^'^      Such   a  representa- 
tion is   given  in  F„    John   [Bj^    p»    lis 

(n+k)/2  c 

5(y)      =      -~ ;r—         \  (y.6j)^  log  y^     d  ij. 


where      log  y  -cj      denotes    the    complex- valued  function 

i^   This  method  was    devised  by  Ro    Courant   and  P.    D,    Lax   [I4.J 
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defined  on  the   real   ajciSo      We    choose      k  =  1      if     n     is 
odd,    and     k  =   0      if     u     is   even,    and   differentiate   ijinder 
the    integral    sign.      This    iiranediately  yields    the   formula— ^i'^ 

(8.12)  5(y)      =      -         "     „       /  log  _  jouJ      do;. 


(2711  j  \uj  1=1 


where   log_   denotes  the  n-th  derivative  of  the  log,  in 
harmony  with  our  previous  notation. 

Let  w(t,y,T,)i^tj)   satisfy  the  following  conditions: 

1.  ^*w  =  0  if   t   <   T, 

2»  '^^'y,l,'C-fhsU)      =  w^(r7y,T,Y,iJ  =  ooo  = 

Ij-o   wCtjV,  :,/..,. j)  =  0  if   t  >  T, 


^  log   ((y-^-|  )  »UJ) 


Prom  the  principle  of  superposition,  we  have 

(80I4)     k(t,y/,n)   =     /      w(t,y,r,V),u;)  du  . 

^  !'-'!=!         ' 

w   can  be  constructed  by  means  of  progesslng  waves, 

provided  that  the  system  J^'^v   =  0   satisfies  the  algebraic 

conditions  listed  in  Section  5o   Moreover,  if  the 

coefficients  of   -C    (and  hence  of  ^^,    are  analytic, 

then  the  formal  series  expansion  for   w   converges  in  a 

neighborhood  of  the  initial  surface  <>   Since  every  solution 

-^    This  formula  was  brought  to  my  attention  by  L.  Girding. 
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of  the  Gauchy  problem  can  be  represented  in  terms  of  k, 
we  are  led  to  a  solution  of  the  Cauchy  problem  by  repeated 
application  of  differentiation^,  integration  and  algebraic 
processeso   We  can  obtain  approximate  solutions  after  a 
finite  number  of  such  operations  by  breaking  off  the 
series  after  a  finite  number  of  termso   This  iSj,  in  a 
sense,  a  generalization  of  the  method  and  results  of 
Hadamard  to  general  analytic  hyperbolic  equationso 

Co   The  generalized  Huygens  principle,, 

The  construction  by  means  of  progressing  waves  shows 
that   wCxy^yCj)   may  be  written  in  the  form 

rak  N_ 
a=T  j=0   J' 


(8.15)  w(x,4,uj)  =  ^^  f.+^,-^i0   U,^,uj))a'^Hx,^,i^)   + 


+  h  (x,C,uJ)  J, 

where   f  (s)  =  log   (s)  „   h   can  be  made  as  smooth  as  we 
o        °=n    -^ 

like  if  N  is  taken  large  enough,  and  the  functions 
^°^(xj4,u))   satisfy  the  characteristic  equation  in  Xo— 2/ 

Lemma  lOolj  to  be  proved  in  Section  10,  implies 
that   k(x<,C)   is  a  smooth  function  of  x   and  £,,      except 
on  the  envelope  with  respect  to  uJ       of  the  surfaces 
^^(x,€,w)  =  0,   ioeo  except  if  x   lies  on  a  ray  which 
passes  through  £,<,      The  locus  of  points  which  can  be 
Joined  to   ?  by  a  ray  is  called  the  ray  conoid  with 

vertex  ^o   The  representation  formiila  (8.8)  immediately 

15/  Note  that  the  characteristic  equations  corresponding 
to  a  system  and  its  adjoint  are  identical. 
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Implies  the 

Generalized  Huygens  Principles'^   The  differenti- 
ability properties  of  the  solution  at  a  point  ^  depend 
on  the  differentiability  properties  of  the  Cauchy  data 
only  at  the  intersection  of  the  ray  conoid  through  4 
with  the  initial  hyperplaneo 

Do   Caustics  and  the  Riemann  functiono 

The  construction  by  means  of  progressing  waves  can  be 
carried  out  only  up  to  the  first  point  where  one  of  the 
phase  siirfaces   0°^{xj,C>(j)  =  0  develops  a  singularity,  i.e. 
the  first  point  of  a  caustico   If  the  principal  part  of 
the  differential  operator  ,^^     has  constant  coefficients, 
then  the  surfaces  jZS  (xj,Cs..J  =  0   are  planes  in  x- space, 
and  no  caustic  can  occur,,   On  the  other  hand^  if  the 
principal  part  of  ^  has  non-constant  coefficients,  then 
in  general  caustics  arise,,   We  could  overcome  this 
difficvilty  by  the  method  of  Section  11,   However,  the 
representation  formula  (808}  also  provides  a  method  for 
extending  solutions  beyond  caustics,,  and  we  shall  use  this 
method  in  order  to  extend  the  Riemann  function  itself 
beyond  caustics. 

Since  the  Riemann  function  is  a  solution  of  the 
differential  equationj,  the  Riemann  function  may  be  repre- 
sented in  terms  of  the  Riemann  function|  it  follows  from 

— /  This  principle  was  first  enunciated  by  Ro  Courant 
and  Po  Do  Lax  [1+]. 
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(8.8)  that 

(8.16)  k(4i,€3)  =  /P°[k(r2,i^,43),  k{^^/r^,Y[)]   di^   . 

This  fact  is  known  as  "Huygens"  major  premise,"—^  or 
the  serai-group  propertyo   Ppovlded  that  Z:"-.      and   T^j 
r^   and   r"^   are  sufficiently  close  to  each  other 
respectively,  the  integrand  in  (8.16)  can  be  constructed 
by  the  method  of  progressing  waves. 

For  the  proof  of  the  generalized  Huygens'  principle 
where  caustics  are  involved,  we  cannot  use  Lemma  10.1 
directlyj  but  we  can  apply  the  same  type  of  reasoning 
to  (8.16).   Prom  (8oli|)  and  (8,15),  it  follows  that 
k(^. ,4*+i)   can  be  expressed  in  the  form 


(8.17)  k(q,q^^)  =  yj^^jjg.z:f'j^.i(i2^%q^p^-))>^ 

^a^^(4^,€j_+iy'-)  +h  i^^s^^^-^suj)  >     duj. 

Equations  (8.16)  and  (8.1?)  show  that  the  smoothness  of 
k(4,,4^)   is  determined  by  the  smoothness  of  integrals 
of  the  form 

(8,18)    J  = 

/'    _   f^^^^(4;j^,C2,^2^^^^^^^%*^2'^2^  ^^2     p^2  ' 

17/  See  J.  Hadamard  L7J,  PP.  53-5^-,  for  a  discussion  of  the 
various  versions  of  Huygens'  principle. 
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where   f  ^^(s)   and  f  ^^(s)   are   C*^  outside  every 
neighborhood  of  the  origin,  and   a     and   a     are 
C°°  functions.  We  denote  the  parameters  in  the  integrand 
by  X  =  (^,, £;-.),   and  denote  the  variables  of  integration 
by  Y  =  i^p'^i^^p"^  °      '^^^^    (80 18)  may  be  abbreviated  to 

(8.19)    J  =  /f^^^j2^^^^X,Y))f^2^jZ^^^^X,Y))a(X,YjdY, 

where      the    integration  is   taken  over   a   certain  manifold 

1  2 

in  Y-space.      Let     S        and     S        be    surfaces    in  Y-space 

given  by  the   equations      0^    ^    =  0      and     0^"^   =  0.      If 

f,      and       TZj,        Z^      and       TT      are    sufficiently   close 

1  2 

to   each  other,      S        and     S        are   regular   surfaceso      This 

is    simply   the   requirement   that   there   be   no   caustics  for 

the  Riemann  functions   in  the    integrand, 

1  2 

(a)    If  S        and     S        are  not    tangent   at    any   of  their 

intersections,    then  we    can  introduce      0  and     0  as 

local    coordinates   in  Y-space    in  the   neighborhoods   of 

1  2 

S        and     S        respectively.      As   in  the  proof   of   Lemma   8.1, 

CO 

we  can  introduce  a  C   partition  of  unity  in  Y-space 

such  that  in  the  support  of  each  function  in  the  partition, 

( 1 )  ( 1 ) 

either  0  r  0  or  0  is  a  local  coordinate,  and 

( 2  )  ( 2  ) 

either  0  "^   Q)      or  0  is  a  local  coordinate.   By 

integration  by  parts  with  respect  to  0  and  0  in 

those  neighborhoods  where  0  and  0  are  local 

coordinates,  we  can  make  the  integrand  as  smooth  as  we  like. 
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(b)    If,    on  the    other  hand,    the   gradients     0y  and 

(2 )  1 

j2L  are   parallel   at   some   point   in  Y-space  where      S 

2  ( 1  )     ( 2 ) 

and  S    intersect,  then  since  obviously  0  =  0  =  0, 

(l )     (2 )  2     "^  1 

we  must  have  0^^      =  0^    '   =   0,   Hence  S,^     lies  on  the  ray 

conoid  through  Cos  and  4^  lies  on  the  ray  conoid 

through  Ct  •   Moreover,  since  0),  is  parallel  to   0i- 

-'•  So  So 

12 
at   the   point   where      S        and     S        Intersect,    these    conoids 

are    tangent  there.      It   follows   from  the   theory   of  first- 
order  partial   differential  equations   that      C-i      and     £,^ 
lie   on  the    same   ray^ 

Hence    the   Riemann  function   is    a   smooth  function  of 
its   arguments   except   if      C,      and      E,^      lie   on  the    same   ray. 
This   is   the    same    as    the   generalized  Huygens  principle. 
A   similar   argument    can  be   used   to   extend  the   generalized 
Huygens   principle   where   an  arbitrary  number   of   intermediate 
points      (lip,    Cp,      etco    are   reqxiiredo 

We  remark   in  passing   that   the    same  reasoning   shows 
that    singularities   In  the    Initial   data  can   only  propagate 
along   characterlsticso      In  fact,    it   would  be   possible   to 
base   the   treatment   of   singularities  upon  the   representa- 
tion formula    (8o8)    instead   of   a  direct   calculation. 
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9.   Preliminary  Remarks  about  Caustics. 

A»   Definition  and  examples  of  caustics. 

We  have  observed  in  Sections  3  and  5  that,  in  general, 
our  formal  construction  cannot  be  carried  out  for  all 
values  of   t,   ¥e  shall  consider  the  case  where  the 
construction  is  impossible  because  one  or  more  of  the 
phase  surfaces  0     =  constant   develop  singularities, 
while  the  system  (5.1)  remains  regular  and  hyperbolic© 
A  connected  set  of  points  where  a  characteristic  surface 
has  such  a  singularity  is  called  a  caustic. 

In  order  to  examine  an  example  of  a  caustic,  let  us 
consider  the  wave  equation  with  n  =  2.   The  characteristic 
surfaces  for  this  equation  can  be  characterized  as  develop- 
able surfaces  whose  rays  are  inclined  at  an  angle  of  l\.S     to 
the  planes   t  =  const.   Thus  the  right  circular  cone  with 
vertex  at  the  point   (1,0,0),   whose  intersection  with  the 
plane   t  =  0   is  the  lonlt  circle,  is  a  characteristic 
surface.   This  siirface  has  a  singularity  at  its  vertex. 
It  can  be  verified  (see  Lemma  9.2  in  this  section),  that 
the  ordinary  differential  equations  (l^..l8)  have  a  singular 
coefficient  at  the  vertex.   Consequently  if,  for  the  wave 
equation,  we  prescribe  Cauchy  data  which  has  a  jump  disconti- 
nuity across  the  unit  circle,  then  the  magnitude  of  the  jiJmp 
in  the  solution  becomes  infinite  at  the  vertex  of  the  cone, 
and  the  formal  construction  breaks  down.   We  may  think  of 
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this  phenomenon  as  a  focusing  of  the  singularity  of  the 
initial  data.   Other  characteristic  surfaces  for  the  wave 
equation  can  be  exhibited  for  which  the  caustic  is  a  curve, 
namely  those  developable s  which  have  an  edge  of  regression. 
In  some  cases,  the  caustic  Itself  may  have  a  cusp. 

Caustics  are  always  involved  If  we  consider  problems 
where  the  initial  data  has  a  "point  singularity"  (for 
example  the  n-dimensional  5-functlon)  j,  since  the  vertex 
of  the  ray  conoid  is  a  caustic  polnt«   If  the  coefficients 
of  the  principal  part  of  the  differential  operator  are 
variable  J,  then  in  general  the  ray  conoid  will  also  have 
singularities  away  from  the  vertex.   This  phenomenon  is 
sometimes  called  "focusing  due  to  the  mediumo"   We  shall 
see  that  there  is  no  essential  difference  between  caustics 
due  to  variable  coefficients  and  caustics  due  to  curvature 
of  the  cxirves  bearing  the  singularity  in  the  initial  data. 

B.   Continuation  of  characteristic  surfaces  beyond  caustics. 

In  spite  of  the  fact  that  the  characteristic  equation 
is  nonlinear,  and  that  caustics  may  have  irregular  shapes, 
the  rays  which  generate  the  characteristic  surface  can 
always  be  continued  beyond  the  caustlco   Hence  the  charac- 
teristic surface  Itself  can  be  continued  beyond  the  caustic, 
This  fact  is  embodied  in 

Lemma  9.1;   In  the  neighborhood  of  a  particular  ray  which 
passes  through  a  caustic,  a  characteristic  surface   S   can 
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be  represented  as  the  envelope  of  regular  characteristic 
surfaces. 

Proof:   Let  I^^   be  the  parameter  along  the  ray  defined 
in  (L(.,l6),  and  let  ^      =■     'C       at  the  caustic  point.   Let 
T^,   be  a  number  less  than   ^  ,   and  sufficiently  close 
to   TT  ,   The  intersection  of  S   with  the  hyperplane 
t  =  t(Tr^)   is  regular.   We  call  this  intersection  "C." 
C   can  be  represented  as  the  envelope  of  its  tangent 
planes o — /   Let   3(lJ  ,  Tl  )   be  the  family  of  character- 
istic sxjrfaces  whose  intersections  with  the  plane   t  =  '^ -. 
are  the  tangent  planes  to   Co   The  theory  of  first  order 
partial  differential  equations  shows  immediately  that  S 
is  the  envelope  of  S(6J,Tr,  )   with  respect  to  6^  .   Since 
the  differential  equation  is  regular  at  ii{'C    ),   the 
surfaces   S(^  ,  ZT^)   are  regular  if   T',   is  sufficiently 
close  to  tT  . 

C,   The  behavior  of  generalized  progressing  wave  solutions 
near  a  caustic. 
We  can  obtain  information  about  the  behavior  of  the 
coefficients   a''   in  the  neighborhood  of  a  caustic  by 
determining  the  nature  of  the  singularity  of  the  coeffi- 
cients of  the  ordinary  differential  equations  (ij-.Sl) 
satisfied  by  the  scalars   cr;^<> 

— '    This  representation  is  given  analytically  in  the  proof 
of  Lemma  10. I4.  in  the  next  section. 
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Lemma  9 o 2 ;   If  a  caustic  of  a  characteristic  surface   S 

intersects  a  ray  which  lies  in  S   at  x(  f  ),   then  along 

this  ray,  the  coefficients  p"^^  are  regular  at  the 
caustic,  and 

(7.1)        P^P(r)  =  pYP  _^_^  ^  pYP(^)  ^ 

where   P'^(')   is  regular  at   L  =  IT  ,   and   p,   is  a 
positive  integer  which  depends  upon  the  shape  of  S   in 
the  neighborhood  of  the  caustico 

Since  this  lemma  is  not  used  in  the  sequel,  we  shall 
omit  the  proofo   The  proof  appeared  in  D.  Ludwig,  Doctoral 
Dissertation,  New  York  University,  December  1958o 

As  an  immediate  consequence  of  this  lemma,  we  obtain 
the  following 

Theorem*   Along  a  ray  which  lies  in  S   and  passes  through 
a  caustic,  the  scalars  cri        can  be  represented  in  the  form 

S-°  =  (r-I^)'^^/^  s  {r)j   where  sA-c)      is  regular  at  r  =  T^. 

Remark;   The  behavior  of  the  coefficients     "^     (j  >  O) 
can  be  determined  in  a  similar  way  from  the  inhomogeneous 
equations  (i4.o21). 

Do   Sijmmary . 

We  can  draw  the  following  conclusions  from  the 
discussion  of  this  sections 

1«   The  fonnal  solution  whose  construction  was 
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described  In  Sections  L\.   and  5  is,  in  general,  not  valid 
for  all  values  of   t,   since  caustics  will  generally  arise 
and  the  coefficients   a-^Cx)   become  infinite  at  caustics. 

2.  Nevertheless,  the  rays  which  generate  the 
characteristic  surfaces  can  be  continued  through  the 
caustic,  and  the  characteristic  surfaces  themselves 

can  be  continued  through  the  caustic  by  taking  envelopes 
of  regular  characteristic  surfaces. 

3.  The  behavior  of  the  coefficients  a"^  near  the 
caustic  can  be  determined  from  the  shape  of  the  charac- 
teristic surfaces  at  the  caustico 
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10 o   Lemmas  on  Integrals  and  Envelopeso 

In  this  section  we  shall  derive  some  properties  of 
integrals  of  distributions  which  depend  upon  parameters. 
The  results  of  this  section  will  make  it  possible  for  us 
to  extend  the  formal  solutions  of  Sections  [4.  and  5  beyond 
caustics,  and  to  investigate  the  behavior  of  solutions  at 
caustics. 

We  can  summarize  the  results  of  this  section  as 
follows s   Let   u(x,4)   be  a  generalized  progressing  wave 
with  the  phase  surfaces   S(4<,c)o   Then  J(x)  =  Ju(x,4)  d£, 
is  a  sum  of  generalized  progressing  waves  whose  phase 
surfaces  are  obtained  by  forming  the  envelope  of   S(4j,c) 
with  respect  to   <|o   The  wave  form  of   J(x)   is  obtained 
from  the  wave  form  of  u(x,4)   by  means  of  fractional 
integrationo   Converselyj,  certain  generalized  progressing 
waves   J(x)   can  be  represented  as  a  superposition  of 
waves  u{x,£,)^    where  the  corresponding  phase  surfaces 
S(£;j,c)   are  planes, 

Ao   Differentiability  properties  of  Integrals  of  waves „ 

As  an  example  of  a  distribution  depending  upon 
parameters,  let  us  consider  the  set  of  functions  Ei0ix^E,)) , 
where   H  is  the  Heaviside  step  function; 


H(s)  =  1    if    s  >  0 
H(s)  =  0    if     s  <  0 
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Here  0{x,£,)      is  a   C°°  function  of  the  vectors  x  and  £,, 
0  (x)  T'  0   if  jZi(x)  =  0.   The  equation  j3(x,4J  =  0   defines 
a  family  of  regular  surfaces   S(4)   in  x-space,  and  the 
function  R{0{x,£,) )      has  a  juinp  discontinuity  across  the 
surface   S(^)o   More  generally  we  may  consider  the  distri- 
bution f(i3J,   where   f(s)   is  a   C°°  function  in  the 
exterior  of  every  neighborhood  of  the  origin.   Now  we  can 
formulate 

Lemma  10.1;  Let  aCx,^)  be  a  C°°  function  of  x  and  ^ 
with  compact  support  in  C~space  for  each  value  of  x,  let 
j3(x5<|)  have  the  properties  described  above,  and  let  J(x) 
be  defined  by 

J(x)  =   f(0(x,^jj  a(x,4)  d£,   , 

where  the  integration  is  taken  over  all  of  4~spaceo   Then 
J(x)   is  a  smooth  function  of   x,   except  along  the 
envelope  of  the  surfaces   S(4)„ 

Proof;   Let  ^Z.   ^^^   be  the  set  of  points  in  4~space 
such  that   0(x,^j  =  Oo   If,  for  a  particular  value  of   x, 
0r  f^  0     for  all   4   in  J[2  (x)  ,   then  0     can  be  intro- 
duced as  a  coordinate  in  ^"Space  in  a  neighborhood  of 
^  (x).   We  can  introduce  a  C   partition  of  unity  in 
4-space  such  that,  in  the  support  of  every  function  in 
the  partition,  either  0   is  a  local  coordinate  or  0  ^  Oo 
The  smoothness  properties  of  J   depend  only  on  those  terms 
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In  the  Integral  where   0   is  a  local  coordinate o   Each  of 
these  terms  may  be  written  In  the  form 


J^  =   J  jf(jZ^)  a^(x,.0,4°)  d0   dc' 


where   a.   has  compact  supporto   By  repeated  Integration 
by  parts  with  respect  to  0^     we  can  make  the  integrand 
have  as  many  derivatives  as  we  likeo   Hence   J  Is 
infinitely  diff erentiable. 

If  J,  on  the  other  hand,  0^   =  0   at  some  point  on 
21  (x),   then  X  lies  on  the  envelope  of  the  siorfaces 
S(C)o   QoEoDo 

Bo   Fractional  integration  and  differentiation  and 
their  adjointSo 
We  can  obtain  more  detailed  information  about  the 
behavior  of   J   in  the  neighborhood  of  the  envelope,  but 
it  is  first  necessary  to  make  some  definitions  and  remarks; 
If  f(s)  is  a  continuous  function  in  the  interval 
-1  <  s  <  Ij   then  we  can  define  integration  of   f  to  a 
fractional  order  by  the  eqiiations 

{10„1)    l"^(a)  f  =  f  (s)   =  — ^^—   [  f(t)  (s-t)"""^  dt 


°-  Pia) 


-1 


(a  >  0  ) 


It  can  be  verified  that   l'*"(a)  l'^(p)  f  =  l"^(a+p)  f, 
¥e  define  the  transformation  T  by  the  equation: 
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(10.2)  T  f(s)   =  f(-s)  o 

We  define  another  transformation  of  f  by  the  equation: 

(10.3)  I~(a)  f  =  T  I°*"(a)  T  f  , 

or 

1 
(10.1+)      l"(a)  f  =  — i—   f  f(t)  (t-s)'^"^  dt  . 

r(a)   J3 

Again,   l"(a)  l"(p)  f  =  l''(a+p)  f  .   It  can  be  verified 
immediately  that  if  f   and  g   are  continuous  then 

(10.5)     /  f^(s)  g(s)  ds   =    (  f(s)  g~(s)  ds  , 
4  ^1 

Thus   I  (a)   and  l'"(a)   are  adjoint  operations.   Fractional 
differentiation  and  its  adjoint  can  be  defined  by  means  of 
a  fractional  integration  followed  by  a  differentiation. 
Fractional  differentiation  and  integration  can  be  defined 
for  distributions  with  support  inside   (-1,1)   by  means 
of  (10.5).   For  distributions  which  are  arbitrarily  smooth 
outside   (-1,1),   fractional  differentiation  and  integra- 
tion may  be  defined  modulo  a  smooth  function  by  means  of 
(I0o5)-   I  (a)   and  I  (a)   do  not  commute,  but  it  can  be 
verified  that   [l'^(|)  l"(|)  -  l'(i)  l'^(~)  J  f   is  a   C°° 
function  if  f   is   C°°  outside  a  closed  subinterval  of 
-1  <  s  <  1  . 
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Co   The  wave  form  and  phase  function  of  the  Integral  of 

a  wave. 

The  following  discussion  is  a  nearly  trivial 

extension  of  the  method  of  stationary  phase. 

Let   f  (s)   be  a  distribution  which  is   C   outside 
o 

a  closed  sublnterval  of   (-Ipl),   Let  the  phase  function 
0{x^£,)      be  a  C°°  function  of  the  vector  x   and  the 
single  parameter  £,0      We  assume  that  0     f-   0.   Hence 
0(x3^)  =  c   Is  the  equation  of  a  family  of  regular  surfaces 
SCc,^)   in  x-spaceo   We  assume  that   ^   is  determined  as  a 
single-valued  function  of  x   by  the  condition  of  station- 
ary phase 

(10.6)  0^(x,C)   =   0  , 

and  write 

(10.7)  €  =  ^(x)o^ 

We  assume  further  that 

(10.8)  <^U,^)      =     P(x,^(x))  +  (4-^(x))P  \|f(x,C)  , 

where  ^{x.^ix))   ^   0.   The  integer   p  may  depend  upon  Xo 
Let   0(x)   be  defined  by 


(10.9)  iZi(x)   s  i2!(x,^(x)) 


=^   If  we  allow  ^(x)   to  be  multivalued;,  slight  modifica- 
tions in  the  following  discussion  are  required. 
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The    STjrface      2(c),      given  by  setting      0{x)    =   c,      is   the 
envelope    of   the    surfaces     S(c,4)»      S(c)      is  regular   if 
i2!^^(x,^(x) )   ^  0    ,      i.e,    if      p  =  2.      Let     a(x,$)      be   a 
C  function  with  compact   support » 

Lemina  10.2:      Let      J(x)      be   defined  by 

(10.10)  J(x)      =  /f^(0(x,C))    a(x,C}    dC    . 

We  consider  three  cases »   Let   +   denote  the  sign 
of   \|/(x,^),   and  +   denote  its  negative. 

(a)  If   p   is  even,  then 

N   — 
(lOoll)        J  =  IZ  f ""  (i2^(yJ)  a^(x)  +  h^(x)  , 

J=0  i+1 

P 
and 

(10.12)  a°(x)   =   +  ^ni/f3)  |^(^^^l"P  a(x,^)  . 

(b)  If   p   is  odd,  then 

N    u.    -^      J..  J 

(10.13)  J  =  1_  [ft,,  (J2((x))  a  ^(x)  +  f"   (0(x))  a"J(x)] 

n-=0    -Hi  Hi 


+  h''(x)  , 

and 

_1 

(10.11+)    a'"°(x)  =  a-°(x)  =  1  r(VP)  I^Cx.rfP  a(x,^)  . 

(c)   If  i^(x,C)   =  ;^(x)  ,   (p  =  oo),   then  it 
follows  at  once  that 

(10.15)  J  =  f  (^(x}]  a(x)  , 

o  ' 
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and 


(10.16) 


'(x)   =   Ja(x,C) 


dC 


In  any  case,  a^(x)  are  C  <,  and  the  remainder 
h  (x)  can  be  made  as  smooth  as  we  like  if  N  is  taken 
large  enough. 

Proofs   It  follows  from  (10.8)  that 

(10.17)   0^(x,€)   =   p(C"^)^°-^  [^(x,C)+4»^)  rego  funct,], 


and 


(10.18) 


J^(T~U 


P  ^(x,^ 


[1+  (4-^)  reg.  functo] 


(10.19) 


Equation  (10.8)  also  implies  that 

^        ,1/P 


<^"^^   =   linhrl  ' '  ^^^  ^^"^^ 


and  hence 


(10.20) 


=      (3gn(g-^)) 


1  =  P 


IKxT^ 


2^ 

!    P 

D"KC-^)reg.functJ 


Now,    introducing     0      as   a  new  variable    of   integration, 
we   have 


(10.21)        J^  = 


(10.22)        J^  = 


il-P. 


a-,(x,4)     \(%ji 


1^ 


^ 


d|2f   , 


iJ2^(x,oo) 


il"P. 


(x,4) 


i2^ 


(sgn(?-^))"°Pf(iZS)   ^^ 


.^1    P 


M 


d^l   p 


where  a^(x,C)   and  &^ix,^)      are  regular  functions.  We 
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note  that 


a. 


Now  we  assTJine  that   p   is  even  and  \|f(x,^)  >  0  ; 
the  proof  for  the  other  cases  Is  obtained  by  obvious 
modification.   Since   \|f(x,^)  >  0   and   p   is  even, 
p(x,co  )   and  ^(x,-oo)  must  both  be  greater  than  0, 
Nothing  is  changed  if  we  set   a,  (Xj,C)  =  0  for 
0   >  0(x,-oo)   and   ap(x,4)  =  0  for  0   >  J2f(x,oo),   Thus 
we  can  extend  the  limits  of  integration  to  the  interval 
^  £  ^  ^  -'••   Combining  J-  and  Jp  , 

(10.23)   J  =  2  f   .f(^)(;^-i^)  ^(^  '"^.  (a,(x,4)  +  a„(x,4)  d0. 

Now  we  expand  the  integrand  in  a  power  series  in  0: 
^^°-^^^  TjTTF  (^1^^^^^  -^  a2lx,4)  -  X^  /^ 

Substituting  (10. 2ij.)  into  (10.23),  we  have 

N        ^ 
(10.25)        J(x)   =  2Z  ^Hl^^^  '^''^''^  ■*"  ^^^""^    » 

which  was  to  be  proved. 

Lemma  10,2  can  be  used  to  prove  the  corresponding 
results  when  4   is  a  vector   ;    4  =  (4  ,..,,4  ).   The 
most  important  case  is  when  j3(x)   is  a  regular  envelope 
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of  0(xj,C),   ioBo  when  p  =  2..   For  this  case,  the  more 
general  result  may  be  stated  as  follows s 
Lemma  10.3'   Let   4  =  (4  jooojC  ),   and  let  f(s), 
0{x,^)    ,       ^(x)   and  a(x, 4)   satisfy  the  hypotheses 

of  Lemma  8.2.   If,  in  addition, 

(10.26)  0  .     .(x,^(x))  <  0  1  =  1,.,.,/, 

(10.27)  0  ^    .(x,^(x))  >  0  i  =  X+l,...,k  , 

then     J(x),      given  by   the    integral 

(10.28)  J(x)      =       j  211  f ^(P(x,4))    a^(x,4)    d4   , 

J    j=o      -^ 

has   the   form   of   a   si:iin  of  waves  plus   a  remainder.      The 
leading   term   in  the   expansion  of      J     is 

(10.29)  Fi0U,KU}}    a(x)    , 
where 

(10.30)  •  F     =     I"    (^)   l'^(|)    f   , 
and 

(10.31)  a(x)      =      (zl/(2.)^/La^(x.^(x))      ^ 

V      i=i  4^4^ 
Proof;   The  lemma  follows  immediately  from  k  times 
application  of  Lemma  10o2. 
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D.   Representation  of  a  family  of  siirfaces  as  the 
envelope  of  a  family  of  planes. 
With  the  aid  of  Lemma  10o3  It  is  easy  to  prove  a 
converse  to  Lemma  10o3s  namely  a  representation  of  certain 
waves  as  a  superposition  of  plane  waveso   First  we  must 
show  how  a  family  of  surfaces  whose  principal  curvatures 
are  nowhere  zero  can  be  represented  as  the  envelope  of  a 
family  of  planes. 

Lemma  10.^.;   Let  x  =  (x  ^ooo^x  )  =  (x  ,,oojX    ,3)  =  (y,s), 
and  let   ^  =  {^^ ,  o  o  o  ,^^~-^  ) .      Let  0(x)   be  a   C°°  fxmction 
of  X   such  that   grad  0(x)  f'  0«   Then  0(x}  =  c   is  the 
equation  of  a  family  of  sxorfaces  in  x-space.   If  none  of 
the  surfaces  S(c)   has  a  zero  principal  curvatiare,  then 
there  exists  a  function  0{x^^)      such  that 

(1)  For  fixed  £,      and  c  ^      the  surface   S(cj,.C)   given 
by  p(x,C)  =  c   is  a  plane. 

(2)  The  planes  S(c,C)   envelope   S(c), 

(3)  ^  ^  -.  7^  0   at  the  envelope   (i  =  lj,«„.^n-l). 
Proof;   Since   grad  0  f   0,   we  can  assume  without  loss 

of  generality  that  0   ^(x  ,<,oo,x  )  7^  0  .   Hence  by  solving 

X 

the   equation 

(10.32)  0(x^,..o,x'^)      =      c    , 

we  can  determine  x^  as  a  function  of   (x  , , . , ,x""  , c ) : 

(10. 33i  x""  =  g(x^,,,.,x''"^,c)  . 
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Equation  (10.33)  can  also  be  written 

(10,314.)  s  =  g(y,c)  . 

Since  we  have  ass-umed  that  none  of  the  principal  curvatures 
of  S(c)   are  zero,  the  matrix 

(g.J  =  (—T-^)  i,j  =  l,,,.,n-l  , 

must  be  non-singular.   By  a  rotation  of  coordinate  axes, 
we  can  insure  that   g..  7^  0   (i  =  l,o«.j,n-l)o 
¥e  define   G(Xj,Cjc)   by  the  equation 

(10035)  G(x,C,c)   =  x^-g{^,c}    -  TZ  sA^,c){x^'£.^)    , 

i=l   ^ 

and  define   S(c,C)   by  setting 

(10036)  G(x,C,c)   =  0„ 

The  surfaces  S(cjC)  are  planes  in  x- space  which  are 
tangent  to  S(c)  at  the  point  {£,  ,.o,j,(g  ,g(^,c))o 
The    equations   for    the    envelope    are 

ir-  n°l  .       . 

(10.37)  ^     s      -  2Z  g^.°^y-  "^    )      =     0        j   =   l,o..,n-l. 

^^-^  i=l      ^^ 

Hence      ^"^   =  x"^      (i   =  lyooo^n-l)      at  the   envelope,    and 

the   equation  of    the   envelope    is 

(10„38)  X      -gCxyooojX        „c)      =     0. 

In  order  to  complete  the  proof  of  the  lemma,  it  is 
only  necessary  to  write  the  equation  for  S(c,4)   in 
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implicit  form,  and  verify  the  condition  for  regularity 
of  the  envelope.   Differentiation  of  (10.35)  shows  that 

(10.39)       ^(x,i,o)  =  .^.g^U^-sl). 

Since     •:)g//')c  =  1/0        ,      ^g/'C   i-   0;      consequently 

n-1  ^^ 
JG/c-)c  7^  0      if     XZ   (x^--^^)^      is    sufficiently   small. 

i=l 
Thus  we    can     determine      c      as    a  function  of      (x,C) 

from   (10o36); 

(10.14.0)  0(x,C)      =      c    . 

A  short  calculation  shows  that,  at  the  envelope, 

(10.1+1)  0      =  -^  j  =  l,.,.,n-l  . 

By  hypothesis,      g..   ^  0    ,      and  hence      <fi    ,    ,   f   Q 
(j   =   l,..<,,n-l).      Q.EoD, 

E,      Representation   of    certain  waves   as   a   superposition  of 

pla.ne   waves. 

Now  we  are  in  a  position  to  show  how  a  distribution 

in  the  form  of  a  wave 

N     -,      . 
(10.14.2)  u(x)   =  J~  F,(?(x))  aJ(x)  , 


J=0  J 

can  be  represented  as  a  superposition  of  plane  waves. 
Lemma  10.5;   Let   P  (s)  ,   0(x)   and  a'^(x)   satisfy  the 
following  conditions; 

(1)  F^(s)   is   C°°  outside  a  subinterval  of   (-1,1) 

(2)  0(x)   is  a   C°°  function,   0   f   0      for  all  x 


x^ 
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such  that   i0(x) I  <  lo 

(3)   a'^(x)   are   C°°  functions  with  compact  support. 

(I4.)   The  principal  curvatures  of  the  surfaces 
j^(x)  =  c   are  different  from  zero  for  all  x   such  that 
|0(x) I  <  1«   /^  of  these  principal  curvatures  have  their 
sign  different  from  the  sign  of  0   o 

Then  u(x)   may  be  represented  as  a  superposition  of 
plane  waves,  plus  remainders.   Specifically,   u  may  be 
represented  as  an  integral  over  (n-1 ) -dimensional  4"Spa-ce 
of  waves  with  the  plane  phase  surfaces  0(x,4)  -   consto 
The  wave  form  of  the  leading  tenn  in  the  integrand  is 
given  by 

Proofs   Let   0(x,4)   be  defined  as  in  the  proof  of 
Lemma  lO.I^o   Let   a°(x,.|)   be  a   C^  function  with 
compact  support  inside  the  region  where   0{x,4)   is 
defined,  such  that 


{»i)V(-i)^  TT  0 


n-l 

^  i  i 
i=l   eg    'io 


(271) 

Lemma  10.3  implies  that 


2 


(x)  . 


(10.1^5)  J'f^(0(x,€))a°(x,C)d€  =  P^(?(x))a°{x)  +  .. 
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where  ...  represents  smoother  terms.  The  next  term 
in  the  expansion  can  be  handled  in  a  similar  way,  and 
so  on. 
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11.   Extension  of  Solutions  Beyond  Caustics. 

¥e  observed  In  Section  9  that  the  generalized  progress- 
ing wave  solutions  described  In  Sections  l\.^    5  and  6  do  not 
In  general  exist  beyond  the  first  point  where  a  caustic  is 
present  ■>   However,  as  we  saw  In  Section  V,  the  family  of 
phase  sixrfaces  of  the  wave  can  be  extended  beyond  caustics 
by  taKlng  envelopes  of  regular  characteristic  surfaces. 
This  fact  is  closely  connected  with  the  theory  of  the 
complete  integral  for  first  order  partial  differential 
equations  (see  R»  Coiorant  and  Do  Hilbert  LlJ,  Chapter  II). 
The  lemmas  of  Section  10  show  that  generalized  progressing 
wave  solutions  have  an  analogous  property;  Integrating  a 
generalized  progressing  wave  with  respect  to  certain 
parameters  £,      corresponds  to  taking  the  envelope  of  the 
phase  surfaces  ^ix,£,)    -   consto   with  respect  to   ^.   Thus, 
with  the  aid  of  Lemma  lOoSj,  we  can  represent  the  solution 
as  a  superposition  of  regular  waves  even  where  caustics  are 
present.   In  the  region  beyond  the  caustic,  the  solution 
can  be  calciilated  as  in  Section  I4.  by  solving  certain 
ordinary  differential  equations  along  the  rays. 

In  order  to  Illustrate  the  method  for  handling 
caustics,  let  us  first  consider  a  simple  example.   We 
assume  that  initial  data  for  the  2-dlmensional  wave 
equation  is  given  which  is  symmetric  about  the  origin. 
Then  the  solution  will  be  a  function  of  t   and  r  (=/x  +y  ) 
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alone.   We  ass\jme  further  that  this  data  is  singular  in 
the  region  1  "^  e  <  r  <  2-e  »   This  singulaj-ity  will 
propagate  along  the  cones  r+t  =  consto,   according  to 
the  theory  of  Sections  2  =  5o   There  is  no  difficulty 
with  the  wave  with  phase  surfaces  r"-t  =  consto,  since 
these  surfaces  are  regular  for   t  >  Oo   On  the  other  hand, 
the  surface  r+t  =  1   is  singxilar  for  t  =  Ij,   and  there 
the  theory  breaks  down.   In  fact,  the  set  of  points  given 
by  r  =  Oj   1  <  t  <  2   is  a  caustic  of  the  family 
r+t  =  c   (1  <  c  <  2)o   (See  Figure  lo ) 

Now,  according  to  Lemma  10.3,  the  initial  data  can  be 
represented  as  a  superposition  of  plane  waveso   The  planes 
in  question  are  tangent  to  the  cones   r+t  =  consto   There 
is  a  corresponding  decomposition  of  the  solution  into 
plane  waves.   Using  the  method  of  stationary  phase 
(Lemma  10,2),  we  can  express  the  solution  as  a  sum  of 
progressing  waveSo   Lemma  10, 2^  also  shows  that,  on  the 
caustic,  the  wave  form  of  the  leading  term  of  the  solution 
is  a   1/2-derivative  of  the  wave  form  of  the  data.   Which 
1/2-derivative  it  is  depends  upon  the  sign  of  0cc(^s^^ 
on  the  cones   r+t  =  const,   before  the  caustic.   After 
passing  through  the  caustic^  the  wave  form  of  the  leading 
term  is  altered  by  the  transformation  I  (3)  I  (-  p")_o   or 
I  (-p)  I  (~  p)  o   Roughly  speaking,  this  a  1/2   differenti- 
ation followed  by  a  1/2  integration,  and  the  differentiability 


91  " 


of  the  solution  is  essentially  unchanged. 

The  following  discussion  is  merely  an  extension  of  the 
method  outlined  above  to  a  more  general  situation. 


Caustic  of 
the  famiTy 


t  =  tT 


\ 


Figure  1 

It  is  sufficient  for  our  purpose  to  be  able  to  carry  out 
the  construction  in  the  neighborhood  of  a  ray  passing 
through  a  point  of  the  caustico   Consider  a  family  of 
characteristic  sxorfaces   S(c)   given  by  setting 
0(x)  =  c,   (c  <  1).   We  assvmie  that  the  whole  family  is 
regular  for   t  =  t^   and  for   t  =  tp,   but  the  caustics 
are  present  in  the  region  given  by   t,  <  t  <  tp* 

We  also  assume  that   tp   is  near   t,,   i.e.  that 
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we  can  represent  the  surfaces   S(c)   as  the  envelopes  of 
a  family  of  characteristic  siirfaces  which  are  regular  in 
the  region  t,  <  t  <  tp  .   An  obvious  way  to  do  this 
would  be  to  represent  the  intersection  of  S(c)   and  the 
surface   t  =  t,   as  an  envelope  of  planes,  and  then 
continue  this  family  of  planes  for   t  >  t,   as  regular 
characteristic  siorfaceso   ¥e  proved  Lemma  10. I4.  with  this 
purpose  in  mind.   Lemma  IO0I4.  states  that  any  regular 
surface  can  be  represented  as  an  envelope  of  planes, 
provided  that  none  of  its  principal  curvatures  are  zero. 
Hov3ver,  the  intersection  of   S(c)   and   t  =  t-,   may  have 
a  zero  principal  curvature »   For  example,  the  converging 
cylinders  \/x  +y  +  t  =  c   in   (tjX,y,z ) -space  have  this 
property.   But  caustics  occur  only  when  certain  principal 
curvatures  of   S(c)   become  infinte.   These  principal 
curvatures  will  certainly  be  different  from  zero  near  the 
caustic.   Thus,  perhaps  using  fewer  than  n-l  parameters, 
we  may  represent  S(c)   as  the  envelope  of  regular 
characteristic  surfaces.   In  analytic  language.  Lemma  lOoi^. 
implies  that  we  can  introduce  parameters   ^  =  (^  ,0.0,4^) 
and  a  function  ^ij,^)      such  that  the  intersection  of  S(c) 
with  the  hyperplane   t  =  t-,   is  the  envelope  of  the 
surfaces  given  by  t(ysi€)  -  c. 

Now  we  can  extend  the  function  t(y,€)   to  functions 
^    (x,4)  =  \|f  (t,y,C)   (ci  =  l,,..,ink),   satisfying  the 
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characteristic  equation  (5o3a)o   This  amounts  to  solving 
the  ordinary  differential  equations  (5o5)  along  the  rays. 
The  theory  of  first-order  partial  differential  equations 
shows  that   S(c)   is  the  envelope  of  certain  of  the 
surfaces   S  (c,C)   given  by  setting  \|/  {x,£,]    =  c.   Moreover, 
t,   and   tp  have  been  chosen  in  such  a  way  that   S  (c,4) 
are  regular  surfaces  in  the  region  t^  <  t  <  tpo 

Assume  that  we  have  a  solution  of  (I^-d)  given  by 

(11.1)  u(x)   =  2ZF"(i^(^))  a^(x)  +  h^(x)  , 

where   u(x)   is  uniformly  smooth  outside  the  region 
|^(x) I  <  lo   In  the  form  (llol),  this  solution  is  only 
valid  in  the  region  below  the  first  caustic  of  the  family 
S(c).   Lemma  10o5  implies  that  we  can  represent  u(t-,y) 
as  a  superposition  of  plane  waveso   The  leading  term  in 
the  integrand  has  the  wave  form 

(11.2)  f   =  I"(^}  I^(-  |)  F^  o 

Here       y(L    is    the    signature    of   the  matrix     \|/    .    ^•(yj^(y)) 
(i,j   <   p)„      Each   of   the  waves    in   the    integrand   can  be 
extended   as   a   solution   of    the    differential    equation  to 
the   region     t-,    <   t  <   tpo 

The   liniqueness    theory   shows   that      u     is    singular   only 
along     S(c)o      We    can   express      u(tp,y)      as   a   sum   of  progress- 
ing waves  by  applying  the  method  of   stationary  phase 
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(Lemma  10o3  )  o 

It  follows  that  the  leading  term  In  the  expansion  of 
u  has  the  wave  form 

where   j^   is  the  signature  of  the  matrix  \j/  .  .(t^jj^^Cy)) 
^ipJ  ^  P ) "   Since  in  general  Ji  ^     X.  s>      ^^^  vave   form  at 
t  =  tp   is  different  from  the  wave  form  at   t  =  t, , 
However,  since   P  is  obtained  from  P   by  p/2   differ- 
entiations  ollowed  by  p/2   integrations,  P  has  essentially 
the  same  differentiability  properties  as   P  „ 

Lemma  10o2  can  be  applied  p   times  to  obtain  an 
expansion  of  u   which  is  valid  at  the  caustiCo   At  the 
caustic,  the  wave  form  will  be  a  certain  fractional  deriva- 
tive of  P  o   At  a  perfect  focus  (as  in  Figure  1)  the  wave 
form  is  given  by   f,,   which  is  obtained  by   p   1/2  differ- 
entiations, and  we  may  say  that   u   loses   p/2   derivatives 
at  the  caustiCo   At  other  types  of  caustics,  less  than 
p/2   derivatives  will  be  losto 
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